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FOREWORD 


This  computer  program  documentation  was  prepared  by  the 
Advanced  Programs  division  of  Rocketdyne,  a  division  of 
North  American  Rockwell  Corporation,  6633  Canoga  Avenue, 
Canoga  Park,  California.  This  document  was  prepared  in 
accordance  with  and  in  partial  fulfillment  of  Contract 
F04611-70-C-0074,  Pulse  Mode  Performance  Model  (Project 
No.  3058,  Program  Element  No.  6.23.02F),  during  the 
period  1  July  1970  to  21  September  1972.  This  contract  was 
administrered  by  the  Air  Force  Rocket  Propulsion  Laboratory 
Edwards,  California.  The  Air  Force  Project  Officer  was 
Capt.  S.  Rosen,  who  replaced  Dr.  Clark  Hawk.  Initially 
Mr.  T.  A.  Coultas  was  the  Rocketdyne  Program  Manager,  with 
Mr,  L.  P.  Combs  replacing  him  just  prior  to  the  program 
extension. 
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INTRODUCTION 


The  Pulse  Mode  Performance  Model  computer  program  has  been  developed 
to  provide  an  analytical  tool  for  accurately  predicting  the  pulse-mode 
performance  of  attitude  control  rocket  engines.  Specifically,  the  prin¬ 
cipal  performance  parameters  predicted  are  propellant  flows,  total  impulse 
and  mean  specific  impulse  for  individual  pulses  and  for  overall  mission 
duty  cycles.  The  pulse  mode  operation  is  applicable  for  pulse  widths 
which  are  long  enough  for  thrust  to  approach  its  steady-state  level  and 
for  pulse  rates  which  are  not  so  rapid  as  to  prevent  thrust  from  decay¬ 
ing  below  10  percent  of  its  steady-state  level  between  pulses. 

This  document,  along  with  the  Fin«.\  Report  (AFRPL-TR-72-16),  contains 
sufficient  descriptive  information  and  instructions  for  knowledgeable 
people  to  use  the  Pulse  Mode  Performance  Model  computer  program  with  a 
minimum  of  difficulty.  This  first  volume  describes  the  computer  program, 
its  required  input  data,  special  operating  instructions  and  output.  Vol¬ 
ume  II  contains  a  listing  of  the  source  program  coding  (excluding  sub¬ 
program  TDK),  of  card  changes  for  special  deck  setups  and  of  the  input 
data  used  in  the  example  case.  Volume  III  contains  the  complete  print 
out  of  the  example  case.  The  last  volume  (IV)  is  a  listing  of  the  TDK 
source  program  coding. 
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\  COMPUTER  PROGRAM  DESCRIPTION 

The  PMPM  compv.ter  program  consists  of  many  subprograms  and  subroutines 
The  complete  "ORTRAN  coding  of  the  PMPM  computer  program  is  presented 
in  Volume  II  and  IV  of  this  program  documentation.  All  of  the  subprograms 
and  subroutines  are  included  in  this  card  listing,  except  for  standard 
mathematical  and  service  routines  which  are  generally  available  in  the 
FORTRAN  library.  In  this  section,  the  primary  purpose  of  each  of  the 
subprograms  and  some  of  the  major  subroutines  are  described  briefly. 

To  aid  in  viewing  the  overall  structure,  execution  order  and  logic,  flow 
charts  of  the  subprograms  (except  for  TDK)  and  some  of  the  major  subrou¬ 
tines  are  presented  in  Figures  1  through  9.  The  complete  description 
of  the  computer  model  is  presented  in  the  PMPM  final  report  (AFRPl-TR- 
72-16). 

MAIN 

The  MAIN  program  of  PMPM  is  an  executive  control  program  which  directs 
the  order  of  executing  its  subroutines  and  subprograms.  Input  control 
data  is  checked  to  determine  which  subprograms  are  to  be  included  in  £he 
analysis.  Subroutine  PPIN  performs  the  function  of  reading  properties 
of  the  propellants  which  are  required  in  one  or  more  subprograms.  Sub¬ 
routine  ENGBAL  solves  the  engine  system  flowrates  and  pressures  based 
on  estimated  performance  efficiencies.  The  primary  subprograms  are  PMDER, 
PULSE  and  DCYCLE. 

PMDER 

The  PMDER  subprogram  block  performs  the  steady-state  combustion  perfor¬ 
mance  analysis.  Its  purpose  is  to  calculate  steady-state  performance 
parameters  such  as:  propellant  flowrates,  chamber  pressure,  characteris¬ 
tic  velocity  (c*),  thrust  coefficient,  thrust,  specific  impulse  and  mean 
fuel  and  oxidiier  combustion  rate  functions.  Subroutine  PMDER  is  the 
executive  control  program  of  this  subprogram  block.  Subprogram:  to  this 
block  include  LISP,  PttSTC,  TRANS  and  TDK. 


LISP 

Subprogram  LISP  performs  the  propellant  injection,  atomize  .ion  and 
spray  distribution  analysis  based  on  specific  injector  .n  para¬ 
meters  and  the  pressure  drop  across  the  injector  orif:-  -s.  Empirical 
spray  atomization  and  distribution  functions  are  bu'it  into  the  pro¬ 
gram  for  common  types  of  impinging  jet  injector  element  types.  Spray 
is  formed  at  the  impingement  points,  and  it  spreads  out  from  each  of 
these  points  in  rays.  The  LISP  analysis  covers  a  region  from  the 
injector  face  to  a  distance  downstream  specified  by  input  data.  At 
the  downstream  location  the  mass  flux  at  each  point  in  a  mesh  plane 
is  calculated  by  summing  the  flux  contribution  from  each  injector 
element.  LISP  defines  the  mass  ar.d  mixture  ratio  spray  distribution, 
estimates  the  amount  of  spray  vaporization  and  predicts  the  mean 
spray  drop  sizes  in  a  plane  downstream  from  the  injector  face. 

PMSTC 

PMSTC  is  the  streamtube  combustion  subprogram  block.  Its  primary 
purpose  is  to  model  the  spray  vaporization  between  the  LISP  region 
and  the  throat  plane.  A  single  stream  tube  analysis  precedes  the 
multiple  stream  tube  analysis  to  obtain  an  approximate  vaporization 
efficiency.  Spray  and  gas  data  from  LISP  are . arranged  into  axisymmetric 
stream  tubes.  Fuel  and  oxidizer  sprays  are  furthor  subdivided  into 
discrete  drop  sizes  to  give  a  distribution  about  their  mean  size. 

The  interaction  of  the  spray  vaporization  and  gas  dynamics  is  calculated 
in  a  step-wise  procedure  downstream  just  through  the  throat.  Stream 
tube  flow  is  coupled  by  the  overall  chamber  area  constraint,  assuming 
constant  pressure  in  axial  planes  until  the  transonic  flow  region 
is  reached.  In  the  later  region,  a  pressure  profile  is  established 
from  a  transonic  flow  solution  from  subprogram  TRANS.  PMSTC  predicts 
combustion  performance  parameters  and  efficiency,  and  it  provide*:  gas 
data  along  a  slightly  supersonic  isobar  tor  the  TDK  nozzle  performance 
subprogram . 
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TRANS 

The  TRANS  subprogram  generates  a  family  of  isobaric  lines  through¬ 
out  the  transonic  flow  regime.  The  isobars  are  used  to  determine 
pressure  profiles  in  the  multiple  stream  tube  analysis  of  PMSTC.  A 
homogeneous  flow  is  assumed  based  on  mean  gas  properties  obtained  from 
the  single  stream  tube  analysis  of  PMSTC. 

TDK 

The  TDK  subprogram  is  a  massive  two-dimensional  (axisymmetric) 
program  which  calculates  the  supersonic,  kinetic  expansion  of  the 
combustion  gas  in  the  nozzle.  Its  purpose  in  PMPM  is  to  provide  the 
thrust  coefficient  efficiency. 

DULSE 

Subprogram  PULSE  characterizes  pulse  performance  by  modeling  the 
transient  performance  of  several  sequences  of  "standard"  width 
pulses  in  which  the  off-time  is  varied  between  pulses  and  by  setting 
up  tables  of  parametric  performance  data.  Its  transient  combustion 
analysis  is  done  in  subroutine  TC0MB. 

TC0MB 

The  TC0MB  subroutine  analytically  simulates  the  transient  performance 
of  a  propellant  feed  system,  propellant  ignition  process,  spray 
vaporization  and  gas  accumulation  in  and  flow  through  the  combustion 
chamber.  Spray  vaporization  is  performed  in  subroutine  GASGEN. 


GASGEN 

t' 

In  subroutine  GASGEN,  fuel  and  oxidizer  spray  ensembles  are  formed 
as  propellant  is  injected  into  the  chamber.  Each  spray  ensemble  is 
vaporized  at  the  rate  specified  by  the  spray  vaporization  rate  functions 
calculated  in  PMS1C.  The  purpose  of  this  subroutine  is  to  determine 
the  transient  flowrates  of  fuel  and  oxidizer  gases  being  supplied  to 
the  chamber. 

DCYCLE 

The  DCYCLE  subprogram  analyzes  the  pul se-by -pulse  and  cumulative 
performance  of  any  sequence  of  pulses,  or  mission  duty  cycle.  Per¬ 
formance  of  individual  pulses  is  synthesized  by  subroutine  SYNTHE. 

SYNTHE 

Subroutine  SYNTHE  synthesizes  pulse  performance  by  using  the  parametric 
performance  tables  of  "standard"  width  pulses  generated  in  PULSE  to 
determine  the  start  and  decay  performance  of  each  pulse  and  by  accounting 
for  pulse  width  by  adjusting  the  center  of  the  pulse  using  steady-state 
performance  for  a  duration  equal  to  the  difference  between  its  width 
and  the  "standard"  pulse  width.  Chamber  wall  temperature  is  calculated 
as  a  function  of  "on-time"  and  "off-time",  and  the  pulse  performance 
is  adjusted  as  a  function  of  wall  temperature. 
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Figure  1.  Flow  Chart  of  PMPM  Control  Program 
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Figure  3.  Flow  Chart  of  LISP  Subprogram 
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Figure  4.  Simplified  Flow  Chart  of  PMSTC  Subprogram 
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figure  5.  Simplified  Flow  Chart  of  PULSE  Subprogram 
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Figure  6.  Simplified  Flow  Chart  of  TC0MS  Subroutine 
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(DECREASE  IN  ENSEMBLE  MASS), 

VOlO  OUT  EHSEmJLES  REDUCED  TQ  MASS 
FRACTION  OF  VAPORIZATION  INEFFICIENCY 


DETERMINE  EARLIEST  time  of 
FUEL/OXIDIZER  IMPINGEMENT. 
WHEN  CURRENT  TIME  REACHES 
IMPINGEMENT  PLUS  IGNITION 
DELAY  TIMES  ANO  MIXTURE  RATIO 
IS  COMBUSTIBLE,  THEN  IGNITE 
GAS 


RETURN 

END 


Figure  ?•  Sinplified  Flow  Chart  of  GASGEN  Subroutine 
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Figure  8.  Flow  Chart  of  DCYCLE  Subprograa 
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Figure  9.  Flow  Chart  of  SYNTME  Subroutine 


DESCRIPTION  OF  INPUT  DATA 


In  this  section,  the  specific  input  data  required  to  execute  the 
computer  program  is  described.  The  input  data  specifications  are  first 
listed,  and  then  some  of  the  more  complex  input  data  is  discussed  in 
detail . 

SPECIFICATION  OF  INPUT  DATA 

The  input  data  entries  required  to  execute  the  computer  program  are 
specified  in  Table  1.  In  this  table,  the  entries  are  listed  by  card 
and  then  further  grouped  by  the  component  model  which  reads  the  data. 
Card  identification  numbers  are  tabulated  in  the  first  column. 

Although  tne  use  of  these  numbers  is  not  mandatory,  card  numbering  per* 
mits  mechanical  sorting  of  the  data  deck,  and  consistency  in  numbering 
aids  in  locating  specific  data  cards.  The  format  for  entering  the 
data  is  also  included  in  tho  first  column,  enclosed  in  parenthesis. 
Standard  F0RTRAN  notatim  is  used  for  the  format  specification. 

The  second  column  of  Table  1  contains  the  FORTRAN  code  names  of  data 
entries.  When  a  parenthesis  enclosing  an  I  and/or  J  follows  the  coded 
name,  then  thi.-  indicates  an  array,  rather  than  a  single  value,  is  re¬ 
quired.  A  description  of  each  entry  appears  along  side  of  each  input 
code  name.  The  description  includes, when  applicable,  such  information 
as:  values  of  indicators  for  program  options,  number  of  cards  required, 
number  of  values  required,  sue  limits  and  dimensional  units. 


Input  data  work  sheets  are  printed  in  Table  2  to  assist  the  user  in 
coding  data  for  key  punching  and  for  providing  a  form  suitable  for 
documenting  specific  input  data  used  on  each  computer  execution  case. 

The  intention  is  for  the  user  to  use  these  worksheets  as  masters  and 
to  reproduce  them  as  they  are  needed. 

The  work  sheets  are  laid  out  for  an  entry  field  width  of  12 
characters  which  is  the  length  generally  specified  in  the  formats.  The 
last  8  character  spaces  on  each  card  are  reserved  for  the  sequence  number, 
and  the  program  does  not  read  these  spaces.  The  FORTRAN  code  name 
appears  to  the  right  of  each  entry  field.  Notes  are  provided  to  in¬ 
struct  the  user  in  situations  in  which  the  inputs  are  to  be  entered  only 
under  certain  conditions.  When  the  note  is  placed  at  the  top  of  the 
work  sheet,  it  applies  to  all  cards  on  that  sheet.  Otherwise,  the 
notes  apply  to  the  specific  card  where  they  appear. 
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TABLE  1.  SPECIFICATION’  OF  INPUT  DATA 


CARD  NO. 

AND  VARIABLE 

FORMAT  CODE 


DESCRIPTION  * 


10,20,30,40  AMAT 
(18A4) 


Identi fication  and  Primary  Control  Data  (HEADER) 
PMPM  case  description.  Four  cards  required. 


50 

(6112) 


60 

(6112) 


ILISP 

ISTC 

ITDK 

IFULSE 

IDCYCL 


JLISP 

JSTC 

JPULSE 

JB0IL 

JIGN 


Primary  control  options  for  selecting  major  subprograms 
LISP,  STC,  TDK,  PULSE  and  DCYCLE,  respectively,  for 
(execution. 

"I"  to  execute,  '‘0"  to  bypass. 

Special:  ISTC  =-l  for  single  stream  tube  (SST) 

'analysis  without  MST  (multiple)  analysis. 

Control  indicators  for  specifying  the  inclusion  of 
propellant  property  data  required  by  subprograms 
(LISP,  STC,  PULSE,  FOIL  and  1GN  respectively. 

"1"  to  include,  "0"  to  omit. 

Special:  JLISP  -  2  required  when  injector  types 
*1,  2  or  3  specified  in  LISP. 

Propellant  Property  Data  (P?IN)  (Cards  90-980) 


90 

( 18A4) 

100 

(6112) 


iu,2,j** 

(6E12.8) 

U0 

(6E12.8) 


211,2,3** 

(6EU.3) 


TITLEP  J  Propellant  description.  One  card. 


NMR  Propellant  property  table  array  sizes  fc-r  mixture  ratio, 

NMACH  mach  number,  nozzle  expansion  area  ratio  and  droplet 

NEPS  film  temperature.  Maximum  sizes:  18,  3,  6  and  20, 

NTK  j respectively. 

Cards  111-573  contain  theoretical  equilibrium  propellant 
combustion  performance  data.  Mixture  ratio  (weight 
oxidizer  to  fuel)  is  the  primary  independent  parameter 
with  mach  number  as  a  secondary  indcpenden:  parameter  on 
chamber  properties. 

THR(t)  Combustion  gas  mixture  ratio  array.  Enter  NMR  values 
arranged  in  either  ascending  or  descending  order. 

TMACH(J)  Combustion  gas  mach  number  array.  Enter  NMACH  values  in 
ascending  order.  Normally,  enter  3  values  with  1st  «  0. 
Last  (or  only)  value  should  be  *  1. 

TSTAT(l)  Combustion  gas  static  temperature  array.  Enter  values  to 
correspond  with  tne  IMR  mixture  ratio  array  and  mach 
number  TMACH(l),  SHR  values  required.  Units:  R. 


*  See  section  on  Discussion  of  Input  Data  for  additional  detail. 

**  Oait  cards  not  required  for  specific  array  size. 
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CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

. DESCRIPTION 

221,2,3* 

(6K12.8) 

TVIS (I , 1) 

Combustion  gas  static  viscosity  array.  Enter  values  to 

correspond  with  the  TMR  mixture  ratio  array  and  madv 
number  TMACH(l).  NMR  values  required.  Units:  lb^/ft-hr. 

231,2,3* 

(6E12.8) 

TGAM(I,1) 

Same  as  above  for  frozen  gamma  array.  Unit  dimensionless. 

* 

241,2,3 

(6E12.8) 

TIM  (1,1) 

Same  as  above  for  molecular  weight  array*  Units:  lb  /mole 

m 

251,2,3* 

(6E12..8) 

TSVON(I) 

Same  as  above  for  sonic  velocity  array.  Units:  ft/sec. 

311-353 

(6E12.8) 

Combustion  gas  property  arrays  corresponding  with  cards 
211-253,  except  for  values  corresponding  with  mach  number 
TMACH(2).  Omit  card  sequence,  if  NMACH  less  than  2. 

411-453 

(6E12.8) 

Same  as  above  except  for  values  corresponding  with  mach 
number  TMACH(3).  Omit  card  sequence  if  NMACH  less  than  3. 

501.,  2, 3* 

( 612 12.8) 

CSTR(I) 

Combustion  gas,  theoretical  c*  array  arranged  to  correspond 
with  TMR  mixture  ratio  array. 

508 

(6E12.8) 

TEPS(J) 

Nozzle  expansion  area  ratio  array.  NEPS  values 
arranged  in  ascending  order. 

511,2,3* 

(6E12.3) 

TCF(I.l) 

Theoretical  thrust  coefficient,  C^,,  array  corresponding 
with  TMR  mixture  ratio  array  and  expanded  to  TEPS(l) 
area  ratio. 

521,2,3 

TOr (1,2) 

Same  as  above  except  for  expansion  to  TEPS(2).  Omit  cards 
if  NEPS  <  2.  . 

531,2,3 

rCF(I,3) 

Same  as  above  except  for  expansion  to  TEPS(3).  Oraft  cards 
if  NEPS  <3. 

541,2,3 

TCF (1,4) 

Same  ns  above,  except  for  expansion  to  TEPS(4).  Omit  cards 
if  NEPS  <4. 

551,2,3 

TCF(I,5) 

Same  as  above  except  for  expansion  to  TEPS(5'>.  Omit  cards 
if  NEPS  <  5. 

561,2,3 

TCF(I,6) 

Same  as  above  except  for  expansion  to  TEPS(6).  Omit  cards 
if  NEPS  <  6. 

*  Oat t  cards  not  required  for  specific  array  size 


Tps^FW  *.  TTSg;-v> ’ 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

570 

(6E12.8) 

DENS1 

DENS  2 

LISP  Propellant  Data  (Cards  570-590.  Include  when 

JLISP^O) 

Density  of  injected  fuel.  Units:  lb  /ft  . 

Density  of  injected  oxidizer-  m 

580 

(6E12.8) 

PX 

DNS AX 1 
DNSAX2 
SDNSA1 
SDNSA2 

CKP1 

Reference  pressure  for  wet  bulb  density.  Units:  psia. 
Saturation  density  of  fuel.  Units:  lb  /ft  . 

Saturation  density  of  oxidizer.  m 

Slope  of  fuel  wet  bulb  density  vs.  pressure. 

Slope,  of  oxidizer  wet  bulb  density  vs.  pressure  * 

Mean  fuel  evaporation  coefficient  (Modified  k1  ) 
in  impingement- fc-collection  plane  region.  f 

Units:  in  /sec. 

590 

(6E12.8) 

. 

CKP2 

STEN1 

STEN2 

VISC1 

VISC2 

Mean  oxidizer  evaporation  coefficient  (modified  V) 
in  impingement- to-collection  plane  region.  .  ° 

Units:  in  /sec. 

Fuel  surface  tension.  Units:  dyne/cm. 

Oxidizer  surface  tension. 

Fuel  viscosity.  Units:  lb  /ft-sec. 

Oxidizer  viscosity.  m 

600.610* 

(6E12.8) 

TVF(I) 

STC  Propellant  Data  (Cards  600-73G.  Include  when  JSTC0O) 

Droplet  film  temperature  array  ranging  from  drop 
temperature  to  maximum  combustion  gas  temperature  used 
in  tabulating  CPVAPF  &  TO/WF,  Enter  NTK  values.  Units:°R 

620,630* 

(6E12.8) 

CPVAPF(I) 

Fuel  vapor  specific  heat  array  corresponding  with  'TVF 

array.  Units:  BTU'lb  -°R. 

a 

540,650* 

(6E12.8) 

TC0NVF(I> 

Thermal  conductivity  of  fuel  vapor- combust ion  gas  film 
corresponding  with  TVF  stray.  Units:  BTU/f t-hr-°R. 

660,670* 

(6E12.8) 

TV0<I) 

Same  as  TVF  except  used  in  tabulating  CPVAP0  &  TCONV0. 

580,690* 

(6E12.8) 

CPVAPjJ(I) 

Same  as  CPVAPF  except  for  oxidizer. 

700,710* 

(6E12.8) 

rcdN*v0(i) 

Same  as  TC0NVF  except  for  oxidizer. 

720 

(6E12.8) 

WBF 

Fuel  normal  boiling  point  temperature.  Units:  °R. 

*  Omits  cards  not  required  for  specific  array  size. 


|,  | 

i  ? 
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CARD  NO. 
AND 
FORMAT 

- SABI 

VARIABLE 

CODE 

p  |  ort.wiriu\i  iun  ur  iixrui  i/run 

DESCRIPTION 

720 

Continued 

TNB/S  ‘ 
RHOI’NB 

Fuel  normal  boiling  point  temperature.  Units:  °R  . 

Fuel  liquid  density  of  at  normal  boiling  point. 

Units:  lb  /ft  . 
m 

Same  as  above  except  for  oxidizer. 

Fuel  critical  temperature.  Units:  R 

Oxidizer  critical  temperature. 

5 

K 

RHOONB 

TCRITF 

TCRIT0 

730 

(6E12.8) 

TBF 

TB0 

RH0BF 

HH0B0 

Fuel  droplet  saturation  temperature  at  P  .  Units:  °R. 
Oxidizer  droplet  saturation  temperature  at  P  . 

Fuel  droplet  density  at  saturation  temperature 
corresponding  to  P  .  Units:  lb  /ft^. 

Same  >>„  RH0BF  except  for  oxidizer. 

740 

(6E12.8) 

WTMLLF 

WTMLI.0 

WTMLVF 

WTMLV0 

DHVF 

DI1V0 

STC,  BOIL  AND  IGN  Propellant  Data (Card  740.  Include 

when  JSTCdO  or  JBOIL^O  or  JIGN^O) 

Fuel  liquid  molecular  weight.  Urits:  lb  /mole. 

Oxidizer  liquid  molecular  weight.  Units:  lb  'mole. 

Fuel  vapor  molecular  weight.  Units:  lb^/mole  . 

Oxidizer  vapor  molecular  weight.  Units:  Ib^mole. 

Fuel  latent  heat  of  vaporization  (plus  sensitive  heat 
to  raise  temperature  of  injected  fuel  to  saturation 
temperature).  Units:  Btu/lb  * 

Same  as  DHVF  except  for  oxidizer. 

/■^Ni 

ON  \Sl 

^.o 

fo 

NRH0F 

NRH00 

NPVAPF 

NPVAP0 

PULSE  Propellant  Datacards  750-980. Include  when  JPUISE^O) 
Size  of  fuel  liquid  density- temperature  table  ( <  20) . 

Size  of  oxidizer  liquid  density- temperature  table  ( <  20) . 
Size  of  fuel  vapor  pressure-temperature  table  (s20). 

Size  of  oxidizer  vapor  pressure- temperature  table  ( s  20) . 

760,770* 

TTRLF(I) 

Temperature  array  for  fuel  liquid  density  array. 

Enter  NRH0F  values.  Units:  R, 

780,790* 

(6K12.8) 

TRU0F ( I) 

Fuel  liquid  density  array  corresponding  with  TTRLF 
temperature  array.  Units:  lb^'ft^ 

800,810* 

(6E12.8) 

TTRLO(I) 

Temperature  array  for  oxidizer  liquid  dansity  array. 

Enter  NRH00  values.  Units:  °R, 

820,830* 

(6KI2.8) 

TRI!00(I) 

Oxidizer  liquid  density  array  corresponding  with 

TTRL0  temperature  array.  Units:  lb  /ft-*, 

OX 

87*0, 8  50* 
(6K12.8) 

trpvF(i) 

Temperature  array  for  fuel  vapor  pressure  array. 

Enter  NPVAPF  values.  Units:  °R. 

*Oralt  or  add  cards  as  required  for  specific  array  size. 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

.  DESCRIPTION 

860,870* 

TPVA?F(I) 

Fuel  vapor  pressure  array  corresponding  with  TTPVF  t • 

(6E12.8) 

temperature  array.  Units:  psia 

880,890* 

TTPV0(I) 

Temperature  array  for'  oxidizer  vapor  pressure  arruy-. 

(6E12.8) 

Enter  NPVAP0  values.  Units:  °R. 

900,910* 

TPVAP0(I) 

Oxidizer  vapor  pressure  array  corresponding  with 

(6E12.8) 

TTPV0  temperature  array.  Units:  psia 

B0IL  and  IGN  Propellant  Data (Card  920.  Include  when 

JB0ILft>  or  JIGN^O) 

920 

TDFPf 

Fuel  freezing  point  temperature.  Units:  R, 

(6E12.8) 

CPLF 

Fuel  liquid  specific  heat.  Units:  Btu/lb  -°R. 

TDFP0 

Oxidizer  freezing  point  temperature.  Units:  °R 

CPL0 

Oxidizer  liquid  specific  heat.  Units:  Btu/lb^-OR. 

IGN  Propellant  Data  (Cards  930-980.  Include  when  JIGN0O' 

930 

CPSF 

Fuel  solid  specific  heat.  Units:  Btu/lb  -°R. 

Fuel  vapor  specific  heat.  Units:  Btu/lb^-OR. 

Fuel  vapor  viscosity.  Units:  lb  /ft-sec. 

(6E12.8) 

CPVF 

MUVF 

KCF 

Fuel  vapor  thermal  Conductivity. 

Units:  Btu/sec-ft^-OR/ft. 

LMBDSF 

Fuel  latent  heat  of  sublimation.  Units:  Btu/lbm. 

LMBDFF 

Fuel  lament  heat  of  fusion.  Units:  Btu/lb  . 

m 

940 

ALPHAF 

Fuel  accommodation  coefficient. 

(6E12.8) 

STENF 

Fuel  surface  tension.  Units:  lb_/ft. 

X 

950 

CPS0 

1  Same  as  on  card  .930  except  for  oxidizer  parameters. 

(6E12.8) 

CPV0 

MUV0 

. 

KC0 

LMBDS0 

LMBDF0 

960 

ALPHA0 

1  Same  as  on  card  940  except  for  oxidizer  parameters. 

(6E12.8) 

STKN0 

f. 

970 

0FINT 

Intermediate  oxidizer/fuel  mass  ratio. 

(6E12.8) 

AIIMI 

Arrhenius  pre-exponential  factor  for  intermediate  times 
tha  molecular  weight  of  the  intermediate. 

Units:  cc-gr/g-mole^-sec « 

EINT 

Activation  energy  for  intermediate.  Units:  cal'g-mgle. 

AQIGN 

Ignition  heat  release  factor.  Units:  Cal-cc/g-tnole  -sec. 

QEXICN 

External  heat  source  for  ignition.  Units:  cal/g-raoic. 

*Uoi»  or  add  cards  as  case  requires. 
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TABLE 

.  SPECIFICATION  OF  INPUT  DATA  (Continued)  • 

CARD  NO. 

AND 

FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

980 

(6E12.8) 

DELHRC 

Heat  of  reaction  for  formation  of  liquid  intermediate. 
Units:  Btu^lb  . 

DELHRV 

Heat  of  reaction  for  formation 

Units:  Btu/lb  . 

m 

of  vapor  intermediate.  i 

TINT 

Temperature  limit  for  formation  of  intermediate. 

Units:  °R . 

TCV 

Temperature  above  which  vapor  intermediate  forms. 

Units:  °R  • 

Engine  Balance  Data  (Cards  1000-1070) 

1000 

NTYPEB 

Type  of  engine  balance i 

(6112) 

"1"  for  constant  pressures  at  feed  system  inlets, 

"2'*  for  constant  injector  end  pressure,  ,  and 

propellant  injection  mixture  ratio. 

NDIA 

Number  of  fuel-oxidizer  orifice  groups. 

1010 

NIF 

Number  of  fuel  orifices  for  this  group. 

(2112, 

NI0 

Number  of  oxidizer  fuel  orificies  for  this  group. 

2E12.8) 

DIF 

Diameter  of  fuel  orifice  for  this  group.  Units:  in. 

Dip 

Diameter  of  oxidizer  orifice  for  this  group.  Units:  in. 

1020* 

Repeat  card  1010  for  each  orifice  group. 

1030 

PVALVF 

Fuel  valve  inlet  pressure. 

(6E12.8) 

PVALV0 

Units:  psta. 

Oxidizer  valve  inlet  pressure. 

>  Value  not  used  if  NTYPEB=2 

XMRI 

Injection  oxidizer/fuel  mass 
flowrate  ratio. 

Not  required  for  NTYPEB=l, 
but  used  as  initial 

PIE 

Injector  end  chamber  pressure, 
Units:  psia. 

estimates  if  non-zero. 

RPC  IN 

Injector  end/nozzle  stagnation  pressure  ratio. 

Not  required,  but  used  as  initial  estimate  if  non-zero. 

1040 

RVAPF 

Estimated  fraction  of  fuel  flowrate  in  vapor  or  gas 

(6E12.8) 

RVAP0 

state  at  chamber  throat.  Set  equal  to  "1."  for  lack 
of  estimate. 

Same  as  above  except  for  oxidizer. 

ECSMIX 

Estimated  mass-weighted  c^h*  mixing  efficiency  factor. 

Set  equal  to  "1."  for  lack  of  estimate. 

ECSENR 

Energy  factor  to  account  for  steady-state  losses  such 

• 

as  heat  losses  to  chamber  wall, 
no  energy  loss. 

Value  of  *’ ?. •  **  is  for 

1050 

DT 

I  Chamber  throat  diameter.  Units:  in. 

(6K12.3) 

RR 

Chamber  radius  ratio  of  throat 
throat  open  rnssn-c. 

wall  profile  to 

*  Omit  or  add  carua  as  case  requires. 
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TABLE  1,  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 

ANP 

VARIABLE 

. DESCRIPTION 

FORMAT 

CODE 

1050 

DCHAM 

Chamber  diameter.  Units:  in. 

Continued 

EBT0L 

Engine  balance  tolerance  on  solution  of  RVAPF,  RVAPO, 

ECSMIX  and  RPCIN  expressed  as  allowable  changes  in 

predictor- to-corrector  values.  Suggested  value:  0.002. 

EBT0L2 

Engine  balance  tolerance  on  injection  flowrates  of 
fuel  and  o:<idi2er  for  LISP  iteration, measured  as  fraction 
change  in  predictor- to-corrector  values.  Suggested, 
value:  0.05  « 

1060 

DLF 

Diameter  of  fuel  feed  system  line.  Units:  in. 

(6E12.8) 

RFLF 

Friction  factor  for  fuel  line  flow  used  in  form; 

144  Rf 

AP  „  *s  — r  Units:  dimensionless. 

1  2  Bc  A2  ' 

(When  AP>>AP^,  R^  may  be  O) 

AMF 

■  Cross-sectional  ar'ea  of  fuel^manifold  at  injector 
orifice  entrance.  Units:  in  . 

RFMF 

Same  as  RFLF  except  for  fuel  manifold  • 

RFIF 

Same  as  RFLF  except  for  fuel  injection  orificies. 

CFIF 

Entrance  loss  coefficient,  K,  for  fuel  orificies  used 
in  form : 

AP  „  K  (1/a2  _  1/a2  units: 

e  2gc P  1  m  Dimensionless. 

Option:  if  negative  value,  then  values-^-  —  /  ref 

where  subscript  FS  refers  to  the  entire  fuel  feeS^systera. 
(Order  of  magnitude  value:  K=2X10~^  X  thrust) 

1070 

. 

Same  as  card  1060  except  for  oxidizer. 

(6E12.8) 

LISP  Subprogram  Input  Data  (Cards  2010-4020.  Include 

when  ILISP/O) 

2010,20 

Two  comment  cards  describing  injector  and  specific 

(18A4) 

LISP  run  conditions. 

Note:  Mass  fluxes  in  an  axial  "collection”  plane  are 
calculated  at  each  (r,0)  mesh  point  of  a  sector.  An 
inner  sector,  which  may  coincide  with  the  complete 
sector,  delineates  wall  and  symmetry  boundaries.  The 
complete  circular  cross-section  is  comprised  cf  cither 
repeating  or  mirror  image  inner  sectors. 

2030 

NEL 

Number  of  injector  elements  included  in  LISP  analysis 

(1216) 

NRML 

(  <  60). 

Number  of  equally  spaced  concentric  arcs  in  (r,6) 
''collection"  plane  mesh  system. 

i 
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TABLE  1,  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

2030 

NTHML 

Number  of  equally  spaced  angular  positions  in  (r ,0) 

Continued 

"collection"  plane  mesh  system.  (NTILMLXNRML  £  400). 

NRWALL 

Number  of  mesh  system  radial  increments  from  center  to 
chamber  wall. (  S  20). 

NTHR 

Index  of  mesh  system  right  (clockwise)  angular 
boundary  of  inner  sector.  (  >  1), 

NTHL 

Index  of  mesh  system  left  (counter-clockwise)  angular 
boundary  of  inner  sector.  (  S  NTHML), 

JSYM 

Type  of  symmetry  indicator: 

"1"  for  mirror  image > 

"2"  for  repeating  image. 

NRBAFR 

Length  in  number  of  radial  increments  of  baffle  on 
right  (clockwise)  boundary. 

• 

NRBAFL 

Same  for  left  (counter-clockwise)  boundary. 

NLSPEC 

Number  of  separate  specifications  for  defining  groups 

of  injector  elements.  (  £  10). 

NCRT 

Number  CPT  constant  radius,  mass  flux  vs.  theta  plots. 

IPUN 

Indicator  for  punched  card  output  (in  PMSTC/STAPE) J 
"0"  for  no  punched  cards. 

"1"  for  punched  cards. 

2040 

IPUNR 

Index  of  right  (clockwise)  boundary  for  punched  card 

(1216) 

output. 

IPUNL 

Same  for  left  (counter-clockwise)  boundary. 

KFCRT 

CRT  contour  plot  indicators  for  fuel,  oxidizer 

K0CRT 

and  total  mass  fluxes  and  reduced  fuel  fraction. 

KTCRT 

respectively.  Enter  number  of  contour  levels  or  "0" 

KFFCRT 

for  no  plot.  With  a  minus  sign,  contour  interval 
will  be  rounded-off.  (  £  35.  Recommend  approx.  12). 

2050 

IRCRT(I) 

Index  of  radius  selected  for  CRT  plot  of  mass  flux 

(1216) 

vs.  theta.  Enter  NCRT  values.  Omit  card  if  NCRT=0. 

2060 

DZ0M 

Chamber  diameter  at  the  "collection"  plane.  Units; in. 

(6E12.8) 

DTHETM 

Theta  (angular)  increment  in  "collection"  plane 
(r,8)  mesh  system.  Units:  degrees. 

THETAR 

Right  (clockwise)  theta  boundary  in  "collection"  plane 
mesh  system.  Units:  degrees. 

Z0M 

Axial  (z)  coordinated  at  the  "collection"  plane. 

Units:  in. 

20M2 

If  non-zero,  an  axial  coordinate  of  a  second  "collection" 
plane.  Units:  in. 

20M3 

If  both  Z0M2  and  7.0M3  non-zero,  an  axial  coordinate  of 
a  third  "collection"  plane.  Units:  in. 

The  collection  plane  used  for  the  PMSTC  analysis  i3 
the  last  non-zero  7.0M, 

r»;rw5H^*:3»  »'  jhv  •■wisr: 
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TABLE  1 

SPECIFICATION  OF  INPUT  DATA  (Continued) 

! 

i 

S? 

l 

CARD  NO. 
AND 
FOMT 

VARIABLE 

CODE 

DESCRIPTION 

6>.'  -  ' 
f 

2070 

CDBAR 

Correlation  coefficient  used  as  a  multiplier  on 

1 

& 

|  ■  ' 

1 

6. 

,f 

£, 

(6E12.8) 

program  calculated  mean  drop  sizes. 

Cards  2010-2080  are  a  set  of  specifications  for  a  group 
of  injector  elements,  and  they  must  be  repeated  NLSPEC 
times. 

I- 

2110 

NTYPE 

Element  type  indicator: 

1 

(1216) 

"1"  for  unlike  doublet., 

"2"  for  like  doublet) 

"3"  fcr  like  doublet  pair,. 

S' 

vf 

k 

"4"  for  triplet) 

g- 

"V  for  pentad  (four-on-one). 

ij’ 

£• 

1 

"8"  lor  general  element  with  correlation 
coefficients  supplied  by  program  user. 

1 

i= 

NPR0P1 

Index  of  propellant  f~om  orifice  1} 

(,1"  for  fuel) 

"2"  for  oxidizer. 

i; 

1 

■ 

| 

NPR0P2 

Same  as  NPR0P1  except  tor  orifice  2. 

IDBAR 

Spray  drop  size  indicator: 

"0"  for  calculating  from  built-in  correlations* 

"l’1  for  user  supplied  values  . 

!:• 

i 

),- 

2120 

DIA1 

Diameter  of  injector  orifice(s)  1.  Units:  in. 

(6E12.8) 

DIA2 

Diameter  of  injector  orifice(s)  2.  Units:  in. 

CDDIA1 

Discharge  coefficient  of  orifice(s)  1. 

£• 

CDDIA2 

Discharge  coefficient  of  orifice(s) 

| 

ZE 

Axial  distance  of  injector  element  impingement  point 
from  injector  plane.  Units:  in. 

1 

x- 

ft 

i‘- 

l 

I 

CANE 

Included  angle,  Y^,  of  orifice  1  and  2  axes.  Units: 
degrees. 

Refer  to  the  text  for  a  full  description  of  the 
orientation  of  an  element's  coordinate  system  expressed 
in  terms  of  rotation  angles  (  a ,  p  and  Y  )  from 
a  reference  coordinate  system  orientation. 

? 

1 

2130 

BETA 

Angle  of  rotation,  p  ,  about  the  y-axis.  Units:  degrees. 

[:■ 

l; 

(6E12.8) 

CA.MMA 

Angle  of  rotation,  Y,  about  the  x  -axis.  Units:  degrees. 

The  following  card, 2140,  is  to  be  included  for  type  3, 

I' 

• 

like  doublet  pair,  elements  only. 

»:' 

L’ 

2140 

CAM FAN 

Included  (cant)  angle  of  spray  fans.  Units:  degrees. 

:.• 

V 

(6M12.8) 

SPFAN 

Y-Spacitu;  between  doublet  pair.  Units:  in. 

£ 

1 

SPKL 

X-Spacing  between  doublet  pair.  Units:  in. 

1 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 


2150 

(6E12.8) 


2160 

(6E12.8) 


VARIABLE 

CODE 


DBAR1 


DBAR2 


. DESCRIPTION 

Card  2150  is  to  be  included  only  if  IBBAR/0  or  NTYPE  >  5. 

Mass  median  spray  droplet  diameter  for  propellant 
from  injector  orifice(s)  1.  Units:  in. 

Same  as  above  except  from  injector  orifice (s)  2. 

The  following  cards,  2050-2080,  are  to  be  included 
for  type  8  injector  elements  only. 

Spray  flux  distribution  correlation  coefficients, 
through  Cg,  for  propellant  from  orifice(s)  1.  Refer 
to  Equation  9  in  -the  PMPM  Final  Report. 


2170 

(6E12.8) 


Same  as  card  2160  except  for  propellant  from  orifice(s)  2. 


3010  LS PEC 

;I12,3E12.8)  RADS 

THETAS 


3020-4000*  LSPEC 
(U2.3K12.81  RADK 

THETAS 

ALFA 


1  Spray  flux  distribution  correlation  coefficient  exponents, 
fa  and  b,  from  propellant  orifice(s)  1.  Refer  to 
equation  9  in  .the  PMPM  Final  Report. 

|Sami  except  for  propellant  from  orifice(s)  2‘. 

Preceding  cards,  2110-2180,  are  to  be  repeated  for  each 
NLSPEC  specification  set.  Sequence  second  digit  of 
card  number  for  each  specification  set,  i.e.,  second 
set  from  2210-2280, 

Data  ort  the  following  card,  3010,  must  be  entered  for 
each  injector  element  in  the  LISP  analysis. 

Index  of  injector  element  specification. (1 S LSPEC  <  NLSPEC) 
Chamber  radial  coordinate,  r_,  of  injector  element 
impingement  point.  Units:  in. 

Chamber  angular  coordinates  0^. ,  of  injector  element 
impingement  point.  Units:  degrees. 

Angle  of  rotation,  or,  about  the  2-axis  used  in  defining 
the  orientation  of  an  injector  element  coordinate 
system  with  a  reference  orientation.  Units:  degrees, 

j  Same  as  card  3010  except  for  next  injector  element. 


NEL  cards  required. 


*  Add  or  omit  cards  as  case  requires. 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 

l 

• 

AND 

VARIABLE 

DESCRIPTION 

FORMAT 

CODE 

The  following  cards,  4010  and  4020,  are  included  only 
if  KFCRT,  K0CRT,  KTCRT  or  KFFCRT  ^  0. 

4010 

WlF 

Minimum  and  maximum  values  for  contour  plots  of  fuel. 

(6E12.8) 

W2F 

oxidizer  and  total  mass  fluxes,  respectively. 

W10 

Options:  W2  *  Wl,  finds  max  &  min  values  from  arrays, 

W20 

WIT 

W2  <  Wl,  finds  max  value  from  array. 

W2T 

.  4020 
(6E12.8) 

W1FF 

W2FF 

J  Same  except  for  reduced  fuel  fraction.  (Range:  Q  to  1.) 
PMSTC  Subprogram  Input  Data  (Cards  5010-5840.  Include 

when  ISTC^O) 

5010 

N0Z0N 

Axisymmetric  stream  tubes  are  formed  from  the  LISP 

(6112) 

NSTPZ 

(r,0)  mass  flux  data. .  After  removing  a  specified 
fraction  of  flowrate  for  the  wall  boundary. stream'  tube, 
the  remaining  flowrate  is  separated  into  N0Z0N  radial 
zones  of  approximately  equal  mass  flow  rate.  NSTPZ 

is  the  number  of  stream  tubes  formed  per  zone. 

Recommended  values:  N0Z0N  =  1,2,  or  3;  NSTPZ  >  6. 

NGT 

Limit:  N0Z0NXNSTPZ  =  18. 

Number  of  spray  drop  size  groups.  Set  =  12  if  ILISP^O. 

NGF 

Number  of  fuel  drop  size  groups,  Set  =6  if  ILISP^O. 

NP 

Number  of  equally  spaced  axial  stations  for  stepwise 
calculations  between  and  including  start  and  throat 
planes.  Maximum  number  of  stations,  including  up  to 

25  downstream  of  the  throat,  is  300. 

NAP 

Number  of  points  used  to  define  chamber  geome  try.  (  <  12).. 

5020 

NSSTI 

Maximum  number  of  complete  passes,  marching  from  start 

(6112) 

NMSTI 

plana  to  throat  plane  in  single  stream  tube  analysis, 
allowed  for  converging  within  a  tolerance  on  the 
solution. 

Maximum  number  of  complete  and  partial  passes  marching 
down  the  chamber,  in  multiple  stream  tube  analysis 
allowed  for  converging  within  a  tolerance  on  the 
solution. 

ICRC 

Number  of  corrector  cycle  passes  at  each  axial  station, 
(normally  =>  1) . 

IPRSST 

Axial- station  print  infcrval  of  single  stream  tube  data. 

XPRMST 

Same  except  for  multiple  scream  tube  data. 

TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

- ... . —  — — ■  - 

DESCRIPTION 

5030,40* 

APROF(I,J) 

Combustion  chamber  wall  profile  data  points  specified 

(6E12.3) 

with  NAP  data  pairs  of  axial  distance  from  ’njector 
and  diameter,  respectively.  Set  first  value  (axial 
distance)  =  0.  Last  data  point  must  be  for  throat. 

Chamber  wall  profile  between  next-to-last  and  last 
points  is  constructed  with  radius  ratio  (RR  input  on 
card  1050)  at  throat  and  a  tangent  line.  Units:  both 
in  inches. 

5050 

DEXIT 

Chamber  nozzle  exit  diameter.  Units:  in. 

(6E12.8) 

PAMB 

Ambient  pressure  at  which  nozzle  flox^  is  discharged. 

Units:  psia. 

ECFVAC 

Rocket  engine  thrust  coefficient  efficiency,  . 

Units:  dimensionless  .  r 

ZSTART 

Axial, z  ,  start  plane  for  stream  tube  combustion 
analysis.  Units:  in. 

ZIMPF 

Mean  axial, z  ,  location  of  fuel  impingement  point(s). 

Units:  in. 

ZIMP0 

Same'  as  preceding  for  oxidizer. 

5060 

crtOl 

Allowable  tolerance  on  the  deviation  of  computed 

(6E12.S) 

contraction  ratio  from  unity  at  chamber  threat. 

* 

The  following  cards,  5100-5120  ,  are  included  only 
when  PMSTC  is  not  preceded  with  LISP.  Data  represents 
stream  tube  conditions  at  ZSTART. 

5100 

NST 

Number  of  stream  tubes.  (  19). 

(6112) 

NASEG 

Number  pie  segments  (sectors)  required  to  represent 
the  complete  cross-section  of  chamber.  Used  as  a 
multiplier  on  areas,  flowrates,  etc. 

* 

5110 

Data  on  this  card  is  for  stream  $ube  l  at  ZSTART. 

(6E12.8) 

AREAl(l) 

Cross-sectional  area.  Units:  in*'. 

GASFL(l) 

Gas  flowrate.  Units:  lb  /sec. 

SMRC(l) 

Gas  mixture  ratio  (oxidTser/ fuel  flowrates). 

SNN(l) 

Number  of  LISP  mesh  points.  Entry  may  be  left  blank. 

SR(1> 

Mean  radius.  Units:  in. 

• 

STH(l) 

Angular  position.  Leave  blank. 

*  Add  or  omit  cards  as  case  requires,. 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued)' 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

. DESCRIPTION 

k 

5120 

Data  on  this  card  is  for  spray  group  1  and  2  in  stream 

(6E12.8) 

tube  1  at  ZSTART . 

GWSPR(l.l) 

Flowrate  of  group  1  propellant  spray.  Units:  lb  /sec. 
Velocity  of  group  1  propellant  spray.  Units:  ft/ sec. 

GVELM(1,1) 

CDIAD1 (1 , 1) 

Diameter  of  group  1  propellant  spray.  Units:  in. 

GWSPR(2,1) 

Same  as  preceding  except  for  group  2  propellant  spray. 

GVELD(2, 1) 

GDIAD1(2,1) 

Data  on  card  5120  is  required  for  each  of  the  NGT 
propellant  spray  groups  in  the  stream  tube. 

Data  on  card  5110  along  with  the  propellant  spray  group 
data  on  card  5120  for  each  of  NGT  groups  must  be 
entered  NST  times. 

Option  in  LISP  allows  punching  of  data  listed  for 
card  5110  and  5120. 

TDK  Subprogram  Input  Data  (Cards  5800-5840;  Include 

when  ITDK^O) 

The  start  line  data  arrays  are  setup  and  punched  on 
cards  in  PMSTC.  If  TDK  is  not  immediately  preceded 
bv  PMSTC,  these  data  must  either  be  input  via  namelist 
(See  "long  form"  option  in  TDK  manual)  or  by  using  the 

punched  cards  with  the  auxiliary  TDK  control  program. 

In  either  case,  the  following  data  must  be  supplied 
in  addition. 

5800,10.20* 

Namelist  Name:  PR0PEL 

(Namelist) 

RSTAR 

Throat  radius.  Units:  in. 

EC 

Chamber  contraction  ratio. 

RWT 

Wall  radius  of  curvature  at  throat/throat  radius. 

Rl 

Wall  radius  of  curvature  at  beginning  of  convergence/ 
throat  radius. 

THETAI 

Wall  convergence  half  angle.  Units:  degrees. 

TUETA 

Wall  angle  at  downstream  end  of  throat  wall  radium. 

Units:  degrees. 

FMOL(N) 

Number  of  moles  of  elemont  N  in  iOOg. 

N  Element  Molecular  Weight 

1  Carbon  12.011 

2  Hydrogen  1.008 

3  Oxygen  16.000 

4  Chlorine  35.457 

5  Fluorine  19.000 

6  Nitrogen  14.008 

l 


*  Add  or  omit  cords  as  case  requires. 
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TABLE  l.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

’ 

■ 

5S00,10,20 

0XM0L(N) 

Number  of  moles  of  element  N  in  lOOg  of  oxidizer* 

Continued 

HF(1) 

Enthalpy  cf  fuel.  Units:  cal/g. 

HF(2) 

Enthalpy  of  oxidizer.  Units:  cal/g. 

5830,40* 

Name list  Name :  DATA , 

Namelist 

IWALL 

Wall  contour  indicator.  Set  =  1.  Other  options  listed 
in  TDK  manual. 

RRT 

Wall  radius  of  curvature  of  downstream  side  of  throat/ 

throat  radius. 

THJ 

Exit  cone  half  angle.  Units:  degrees. 

EPS 

Expansion  area  ratio. 

Nl 

Print  interval  of  calculations  for  interior  points 

along  characteristic  lines  selected  for  print. 

N2 

Print  interval  of  calculations  for  characteristic  lines. 

IMAX 

Limit  on  iterations  for  convergence.  15  recommended. 

IMAXF 

Termination  indicator  if  IMAX  not  sufficient. 

"1"  to  terminate.  Recommend  "O"  to  continue. 

PULSE  Subprogram  Input  Data  (Cards  6010-6630.  Include 

when  IPULSE^O) 

6010 

NS PACE 

Number  of  pulses  in  sequence  of  "standard"  width  pulses 

(6112) 

with  spacing  (off-time)  as  a  parameter  (  <  12). 

NTWALL 

Number  of  sequences  of  "standard"  width  pulses  with 
chamber  wall  temperature  effect  on  boil-off  as  a 
parameter.  (  <  6) . 

XIGN 

Ignition  model  Indicator:  "1"  to  envoke,  "0"  to  bypass. 

IBOIL 

Same  except  for  boil-off  model. 

IPPRT 

Print  indicator  for  which  each  level  adds  to  printout: 

"O"  for  no  transient  printout, 

"1"  for  regular  transient  (2  lines  per  At)printout, 

"2"  for  feed  system  data  from  FLjjpW  and  BpIL, 

"3"  for  dynamic  flow  solution  data  from  FL0W. 

Normal  value  is  l, with  2  and  3  used  for  checkout 
only. • 

ICRTP 

CRT  plot  indicator:  "V  to  plot  thrust  and  flowrate 
traces,  "0"  to  omit. 

6020 

DIMS 

Time  step  for  transient  calculations.  Units:  msec. 

(6E12.8) 

STDPW 

Pulse  width  (on- time)  for  "standard"  pulses.  Should 
be  long  enough  for  transients  to  decay  and  as  short 
as  possible  to  minimize  computer  execution  time. 

• 

Units:  msec. 

TAMB 

Ambient  temperature  for  initial  chamber  temperature. 

Units:  °R. 

PAMS 

Ambient  pressure  at  which  nozzle  flow  discharges. 

Units:  psla. 

T*  -« 

*  •*  -1  • 

Tcmerature  of  oxidizer  cnrcrln?  feed  svstem.  Units:  °R. 

*  Add  or  omit  cords  as  ease  requires. 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

• 

DESCRIPTION 

6020 

Continued 

TLF 

Temperature  of  fuel  entering  feed  system.  Units:  °R. 

6030 

VVALF 

Fuel  valve  applied  voltage.  Units:  volts. 

(6E12.8) 

VVAL0 

Oxidizer  valve  applied  voltage.  Units:  Volts, 

C0EHTH 

Exponential  coefficient  for  chamber  wall  heating  during 
continuous  firing. 

C0EHTC 

Exponential  coefficient  for  cooling  chamber  walls 

while  chamber  is  off. 

TAUIGC 

Ignition  delay  constant,  r  .  Units:  msec. 

Steady-state  thrust  charaber°wall  temperature  Units:  R. 

TCWSS 

6040 

QSBFSb 

Steady- state  heat  soak-back  from  fuel  feed  system  to 

(6F.12.8) 

propellant.  Units:  Btu/sec. 

QSB0SS 

Same  except  for  oxidizer. 

TFSBF 

Steady-state  oxidizer  feed  system  temperature. 'Units:  °R 

TFSB0 

Same  except  for  fuel. 

6050,60 

PSPACE(I) 

Pulse  spacing  (off-tine)  array.  Enter  N5PACE  values. 

(6E12.8) 

Units:  msec. 

6070 

TTVTALL(J) 

Chamber  wall  temperature  array.  Enter  NTWALL  values 

(6E12.8) 

ranging  from  ambient  to  steady-state  values.  Units:  °R. 

6100 

NTVEF 

Array  size  for  fuel  valve  opening  table  of  fraction 

(6112) 

open  area  vs.  elapsed  time  from  valve  coil  de-enefgizatior 
(S9). 

NTVDEF 

Array  size  for  fuel  valve  closing  table  of  fraction  open 
area  vs.  elapsed  time  from  value  coil  de-energization 
(  5  9). 

NTVEO 

Same  as  NTVEF  except  for  oxidizer  valve. 

NTVDEO 

Same  as  NTVDEF  except  for  oxidizer  valve. 

it 

6110 

TTVEF(I) 

Time  array  for  fuel  valve  opening  table.  Enter  NTVEF 

(6E12.3) 

values.  Units:  msec. 

6120* 

TAVEF(I) 

Fraction  open  area  array  for  fuel  valve  opening  table. 

(6E12.8) 

Enter  NTVEF  values. 

6130* 

TTVDEF(i) 

Time  array  for  fuel  valve  closing  table.  Enter  NTVDEF 

(6E12.8) 

values.  Units:  msec. 

6140* 

TAVDEF(I) 

Fraction  open  area  array  for  fuel  valve  closing  table. 

(6E12.8) 

Enter  NfVOEF  values. 

6150* 

TTVK0(I) 

Same  as  TTVEF  except  for  oxidizer  table. 

(6E12.3) 

*  Add  or  omit  curds  as  case  requires. 
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CARD  NO. 
AND 
FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

6160* 

(6E12.8) 

TAVE0(I) 

Same  as  TAVEF  except  for  oxidizer  table. 

k 

6170 

(6E12.8) 

TTVDE0(I) 

Same  as  TTVDEF  except  for  oxidizer  table. 

6180* 

(6E12.8) 

TAVDE0(I) 

Same  as  TAVDEF  except  for  oxidizer  table. 

6190 

(6E12.8) 

AVEF 

CFVF 

2 

Fuel  valve  full  open  area.  Units:  in  • 

Entrance  loss  coefficient,  K,  for  partially  open 
fuel  value  used  in  form: 

. _  144  w2  ..  /  1  1  \ 

4P  VV  '  A2  ). 

'  open  7 

Units:  dimensionless. 

AOVEF 

A1VEF 

AOVDEF 

AlVDfTF 

)  Coefficients  and  a^  for  fuel  valve  energize 

)  time  calculation;  t  =  a  +  aiv.  Units  of  t:  msec. 

}  Same  as  preceding  except  for  valve  de-eneegize  time. 

6200 

(6E12.8) 

AVE0 

CFV0 

AOVE0 

A1VE0 

AOVDE0 

A1VDE0 

Same  as  card  6190  except  for  oxidizer  valve. 

6210 

(6E12.8) 

ILF 

LL0 

V0LMF 

V0LM0 

LIF 

LI0 

Fuel  feed  system  line  length.  Units:  in,  (>0.) 

Oxidizer  feed  system  line  lengthy  (>0.) 

Fuel  manifold  volume.  Units:  in  .  (>Q.) 

Oxidizer  manifold  volume.  (>0.) 

Moan  length  of  fuel  injector  orifice(s>.  Units:  in.  (>0.) 
Mean  length  of  oxidizer  injector  orifice(s).  (>0.) 

Cards  6300-6500  are  included  only  if  PULSE  is  not 
preceded  by  PMSTC  in  the  same  computer  run.CISTC^O) 

6300 

(6E12.3) 

V0LC 

LCHAM 

AEXIT 

AUCHAN 

DIMPF 

DlMI’0 

3 

Combustion  chamber  volume.  Units:  in  « 

Chamber  length.  Units:  in.  0 

Chamber  nozzle  exit  area:  in'".  2 

Chamber  walL  area.  Used  in  IGN  only.  Units:  in  . 

Fuel  injection  impingement  distance  from  injector. 

Units:  in. 

Oxidizer  injection  impingement  distance  from  injcctot. 
in. 

*  Add  or  omit  cards  as  case  requires 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO, 
AND 
FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

6310 

F.CFVAC 

Thrust  coefficient  efficiency  in  vacuum. 

(6E12.8) 

WCCHAM 

Weight  of  gas  in  chamber  during  steady-state  operation. 

Units:  lbm. 

6321-6329 

TR(I) 

Time  array  for  the  depletion  of  propellant  spray 

by  10 

ensembles.  Enter  50  equally  incremented  times. 

(6K12.8) 

These  data  are  punched  when  PMSTC  is  executed. 

Units:  msec. 

6350-6338 

SPMF(I) 

Depletion  (resulting  from  gasification)  array  for  a 

by  10 

fuel  spray  ensemble  corresponding  with  the  TR  time  . 

(6E12.8) 

array,  expressed  as  a  fraction  of  original  spray 

ensemble  mass  as  a  function  of  elapsed  time  from  the 

instant  of  impingement. 

6339-6347 

SPRAYO(I) 

Same  as  SPRAYF  except  for  oxidizer  ensemble. 

by  10 

(6E12.8) 

Include  cards  6600-6630  only  if  IGN  is  executed. (IIGN^O) 

6600 

LIWF 

Mean  path  length  of  fuel  drop  from  inector  to  chamber 

(6E12.8) 

wall.  Units:  ft. 

LITF 

Mean  path  length  of  fuel  drop  from  injector  to  thoat. 

Units:  ft. 

LIW0 

Same  as  LIWF  except,  for  oxidizer. 

LITp 

Same  as  LITd  except  for  oxidizer. 

;  6610 

LF 

Number  of  passes  for  looping  on  fuel  and  oxidizer 

(6112) 

L0 

sections  of  program.  Set  all  of  them  =  1. 

LIF 

LI0 

6620 

TMAX 

Maximum  simulated  time  allowed  Co  attain  ignition.* 

(6E12.8) 

Units:  msec  (8,)jf 

DPRT 

Increment  on  simulated  time  for  pringtng  output. 

Units:  msec.  (8.)* 

THTF 

Time  interval  during  which  fuel  drop-chamber  wall,  heat 

! 

transfer  factor  changes  from  maximum  to  minimum  value. 

* 

* 

Units:  sec.  (0.002)*. 

j 

FMAXF 

Maximum  and  minimum  fuel  drop/chamber  wall  heat  transfer 

; 

FMINF 

)  factor,  respectively.  Units:  dimensionless. 

f 

/  (0.125  &  0.01)*. 

TUT0 

Same  as  THTF  except  for  oxidizer. 

.  6630 

FMAX0 

Same  as  FMANF  &  FMINF,  respectively,  except  for 

j  (6E12.8) 

FMIN0 

oxidizer.  (0.25  &  0.02)*, 

t 

1 _ 

OTNMS 

IGN  calculation  time  interval.  Units:  msec.  J 

*  Values  used  in  Seamans'  sample  case. 
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TABLE  1.  SPECIFICATION  OF  INPUT  DATA  (Continued) 


CARD  NO. 

AND 

VARIABLE 

. DESCRIPTION 

FORMAT 

CODE 

DCYCLE  Input  Data  (Cards  7400-7440.  Include  when 

IDCYCL^O) 

If  DCYoLK  is  not  preceded  in  same  computer  run  by 

PULSE,  data  cards  punched  by  PULSE  must  be  inserted  here. 

7400,10 

Two  comment  cards  describing  duty  cycle  case. 

(18A4) 

-  '  - 

7420 

NSEQ 

Number  of  pulse  sequences  in  duty  cycle. 

(IU, 

ECFQ 

Thrust  correlation  coefficient  used  in  calculating  a  total 

5E12.8) 

impulse  adjustment  variable  which  is  dependent  on 
chamber  wall  temperature. 

7430 

NPS 

Number  of  pulses  in  first  pulse  sequence. 

(2112, 

JPAGE 

Pri.it  control  indicator  for  full  page  output. 

4E12.8) 

'O’*  to  suppress, 

"l"  for  first  pulse  only, 

"2"  for  first  and  last  pulse, 

"3"  for  first,  center  and  last  pulse, 

"4"  for  all  pul: 2s . 

"5"  Same  as  ”3",  but  also  deletes  all  short  format 
•  printout . 

PWIDTH 

Pulse  width  (duration).  Units:  msec. 

0FFB 

"Off-time"  (spacing)  between  pulses.  Units:  msec. 

0FFC 

"Off- time"  between  last  pulse  and  next  sequence 

For  last  sequence,  enter  maximum  time  in  PSPACE  array. 

Units:  msec. 

7440* 

NPS 

Same  as  card  7430  except  for  2nd  pulse  sequence. 

JPAGE 

PWIDTH 

0FrB 

Provide  separate  card  for  each  of  NSEQ  pulse  sequences. 

END  OF  PM?M  INPUT  DATA 

*  Add  or  omit  cards  as  case  requires. 
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Include  card  if  NMR>12 


or  JPUtSE/0 


JSTC/O  or  JPULSE^O 


JSTCj^O  or  JPULSEj^O 


J$TCj*0  or  JPULSE^O 


(Ms-indr  jo  o^oisr 


ese  cards  If  JSTC/O  or  JPULSf^O 


Include  these  cards  If  JLISP^O  or  JSTCj^O  «  r  JPULSE^O 


these  cards  ff  JSTC^O 


LM8DS0 

LMBDF0 


FOR  NTYPE=3,  ONLY 


N™EF  1  LU  I .  1  TAVDEF  ( I 


JPUISE 


DISCUSSION  OF  INPUT  DATA 


This  section  supplements  the  description  of  input  data  in  Table  1  .  In  general, 
only  information  relative  to  those  entries  which  need  further  explanation  for 
proper  usage  are  discussed  here.  The  discussion  is  grouped  according  to  the 
subprograms  which  read  the  input  data. 

PMPMID  Input  Data 

The  purpose  of  this  group  of  input  data  is  to  identify  the  computer  run  and 
to  specify  primary  control  indicators. 

Computer  Run  Comment  Cards.  Four  comment  cards  are  read  and  printed  for  the 
purpose  of  identifying  the  computer  run  output  for  future  reference.  The 
comments  should  identify  the  rocket  engine  being  simulated,  nominal  operating 
conditions  and  purpose  of  the  analysis.  The  comment  cards  are  printed  on  a 
title  page  which  identifies  the  computer  model,  including  a  revision  date. 

The  date  of  the  computer  run  is  automatically  printed  (i.e.,  if  the  computer 
hardware  is  equipped  to  do  so) . 

l 

Primary  Model  Control  Indicators.  Control  indicators  provide  flexibility  in 
executing  PMPM  as  one  complete  comprehensive  run  or  individually  by  major  sub¬ 
program  blocks.  Allowable  combinations  for  running  subprogram  blocks  are: 


Model 


Indicator 

1 

,■> 

C 

3 

4 

j> 

6 

•**• 

ILISP 

1  or  0 

1 

or  0 

0 

1 

0 

0 

ISTC 

1 

1 

0 

0 

±1 

0 

ITDK 

1  or  0 

1 

or  0 

0 

0 

0 

0 

IPULSE 

1 

0 

1 

0 

0 

0 

IDCYCL 

1  or  0 

0 

1  or  0 

0 

0 

1 

These  control  indicators  are  provided  to  allow  flexibility  in  executing  PMPM. 
Types  of  situations  for  which  they  are  useful  include: 

1.  Execution  of  individual  subprograms  when  the  entire  model  is 
not  required  for  a  particular  analysis 
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dfcriir  ,ic 


2.  Bypass  subprograms  which  cannot  accommodate  a  particular  case 


3.  Execute  additional  cases  without  rerunning  preceeding  sub¬ 
programs 

4.  Enable  the  user  to  evaluate  results  before  proceeding  into 
the  next  subprogram  analysis 

Whereas  PMPM  is  fully  integrated  to  execute  a  complete  analysis  in  a  single 
computer  r?m,  execution  of  the  subprograms  in  sequential  runs  is  generally 
recommended  to  enable  the  user  to  review  the  results  from  one  before  pro¬ 
ceeding  to  the  next.  In  particular,  steady-state  performance  predicted  by 
LISP  and  FMST'C  should  be  inspected  before  proceeding  with  an  expensive  TDK 
run  or  starting  the  pulse-mode  analysis. 

If  TDK  is  included  in  the  analysis,  a  separate  run  is  recommended  using  the 
TDK  auxiliary  control  program  and  the  punched  data  cards  generated  by  PMSTC. 
Execution  of  TDK  from  PMPM  without  preceding  it  with  PMSTC  in  the  same  run, 
is  undesirable  in  that  the  start  line  data  arrays  nua.  be  manually  punched 
using  a  namelirt  format. 

Control  Indicators  for  Propellant  Properties.  In  order  to  limit  the  input 
data  requirements  to  those  required  by  the  subprograms  selected  for  execution, 
indicators  are  read  to  specify  subprograms  for  which  propellant  data  is 
supplied.  Once  a  propellant  data  deck  is  prepared,  it  may  be  convenient  to 
keep  it  intact  even  though  some  parts  cf  it  are  not  required  for  a  particular 
case.  The  indicators,  being  independent  of  the  model  selection  indicators, 
allow  for  this  option.  Checks  are  made  to  terminate  the  case  if  insufficient 
propellant  data  is  requested  relative  to  the  subprogram  selection  indicators. 

PPIN  Input  Data 

Propellant  property  data  required  by  the  different  subprograms  overlap  con¬ 
siderably.  To  eliminate  redundancy,  all  propellant  duta  is  read  by  a  siugle 
subroutine,  PPI\\  before  executing  any  of  the  major  subprograms. 

i 

Combustion  Gas  Properties.  Combustiou  gas  properties,  which  must  be  obtained 
from  prior,  peripheral  computation  using  a  thermodynamic  equilibrium  perfor¬ 
mance  computer  program,  are  entered  in  tabular  form  as  part  of  the  input. 
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Rocketdyne's  free-energy-minimization  program  has  been  used  here  to  determine 
these  properties,  but  any  comparable  program  is  adequate.  The  combustion  gas 
property  table  is  identified  in  the  FMPM  printout  with  a  description  supplied 
by  the  program  user  on  a  single  input  data  card.  This  description  should  spec¬ 
ify  the  propellants  used,  the  pressure  at  which  the  combustion  analysis  was 
made  and  any  special  information  about  the  combustion  analysis.  Static,  shift¬ 
ing  equilibrium  properties  of  temperature,  viscosity,  molecular  weight  and  sonic 
velocity,  and  frozen  specific  heat  ratio,  all  as  functions  of  both  mixture  ratio 
and  Mach  number,  are  required  in  the  table.  Also  theoretical  shifting  perfor¬ 
mance  of  characteristic  velocity  (c*)  as  a  function  of  mixture  ratio  and  of  thrust 
coefficient  (C^)  as  a  function  of  both  mixture  ratio  and  nozzle  expansion  area 
ratio  are  required.  The  mixture  ratio  array  must  cover  nearly  an  unlimited 
range  to  sufficiently  cover  typical  situations  encountered  in  both  LISP  and 
the  multiple  stream  tube  section  of  PMSTC.  In  order  to  cover  the  range  from 
0  to  infinity  and  to  improve  the  interpolation  accuracy,  a  reduced  oxidizer 
fraction  array,  bounded  by  0  and  1,  is  constructed  internally  from  the  mixture 
ratio  array,  c^  and  used  thereafter  in  its  place.  The  normalicing  factor 
is  the  central  value,  c^,  from  the  mixture  ratio  array,  and  is  used  as  follows: 

F*  Ci 

o,i  c.  +  c 
’  x  m 

In  selecting  mixture  ratio  values  for  the  table,  a  plot  of  static  temperature 
vs  F^,  with  near  the  stoichiometric  mixture  ratio,  is  heipful  in  selecting 
values  best  suited  for  linear  interpolation. 

The  Mach  number  array  consists  of  up  to  three  values.  Normally,  the  values 
should  be  0,  a  value  approximately  equal  to  the  Mach  number  at  the  beginning 
of  nozzle  convergence  and  1.  If  limited  to  two  values,  0  and  1  are  recom¬ 
mended.  For  a  single  value,  a  value  of  1  is  required  to  obtain  a  nozzle 
flowrate  corresponding  with  shifting  equilibrium  propellant  performance. 

LISP  Saturation  Density.  Wet  bulb  density  is  required  in  LISP  for  adjusting 
spray  drop  size  to  account  for  a  change  in  propellaut  density  from  injection 
to  vet  bulb  conditions.  Since  wet  bulb  conditions  are  not  known  for  input, 
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saturation  conditions,  which  are  at  a  temperature  usually  only  slightly  higher 
than  wet  bulb,  are  used  instead.  Propellant  property  input  data  required  to 
calculate  saturation  density  of  fuel  and  oxidizer  for  specific  chamber  condi¬ 
tions  are  saturation  densities,  DNSAX1  and  DNSAX2,  at  a  specified  reference 
pressure,  PX,  along  with  their  respective  slopes,  SDNSA1  and  SDNSA2,  of 
saturation  density  with  respect  to  pressure.  .  Pressure  has  little  direct  effect 
on  propellant  density  at  constant  temperature,  but  has  a  significant  effect 
on  saturation  temperature  and  the  density  corresponding  with  it.  The  latter 
effect  is  the  one  which  must  be  accounted  for. 


LISP  Evaporation  Coefficients.  Mean  evaporation  coefficients  CKPl  and  CKP2 

for  fuel  and  oxidizer,  respectively,  are  used  in  approximating  the  partial 

propellant  evaporation  in  the  LISP  region  of  the  combustion  chamber.  Only 

approximate  values  can  be  supplied  for  these  coefficients,  and  values  on  the 
-4 

order  of  4  x  10  may  suffice  in  general.  However,  for  more  specific  values, 
the  relation  should  be  considered  of  these  coefficients  with  specific  values 
of  the  evaporation  coefficients  calculated  in  and  printed  from  subprogram 
PMSTC  on  previous  computer  runs.  Values  of  the  evaporation  coefficients  in¬ 
put  for  use  in  the  LISP  region  generally  should  be  about  l/4  or  l/ 5  of  those 
calculated  for  use  in  the  combustion  region  in  order  to  account  for,  in  parti¬ 
cular,  incomplete  atomization  over  the  LISP  region. 


PMSTC  Propellant  Vapor  Properties .  Tables  of  fuel  and  oxidizer  vapor  speci¬ 
fic  heat  and  thermal  conductivity  as  functions  of  temperature  are  required, 
spanning  the  range  of  temperatures  across  the  vapor/combustion  gas  films 
around  spray  droplets.  At  the  lover  temperatures ,  specific  heats  at  constant 
pressure  may  be  conveniently  obtained  directly  from  propellant  enthalpy  tables 
or  charts.  At  higher  temperatures,  dissociation  is  important  and  it  is  appro¬ 
priate  to  bleud  the  low  temperature,  undissociated  data  into  equilibrium  dis¬ 
sociation  data. 


The  thermal  conductivity  needed  is  uot  simply  that  of  the  vapor,  but  that  of 
the  combustion  gas-vapor  mixture  between  a  droplet's  surface  and  a  surround¬ 
ing  flame-front.  Again,  a  blending  between  undissociated  propellant,  disso¬ 
ciated  propellant  and  propellant-rich  combustion  gases  is  appropriate  at  low 
temperatures  with  a  gradual  shift  to  the  conductivity  of  the  combustion  ga.«.es 
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aloue  at  high  temperatures.  The  general  technique  used  for  generating  values 
for  these  tables  ^Ls  to  plot  thermal  conductivity,  calculated  from  a  thermody¬ 
namic  combustion  perforaiance  computer  model,  as  a  function  of  combustion  gas 
temperature.  Generally,  this  plots  in  the  fe-'m  of  p.  loop  in  vhich  the  upper 
branch  corresponds  with  fu>l-vich  mixture  ratios  and  the  lower  branch  corres¬ 
ponds  with  oxidizer-rich  mixture  ratios.  On  the  same  grid,  thermal  conducti¬ 
vity  of  the  mire  va-poi*  are  also  plotted.  Then  the  conductivity  line  for 
the  film  around  a  fuel  drop  is  constructed  by  starting  along  the  fuel  vapor 
line  at  low  temperatures  and  blending  it  with  the  fuel-rich  branch  of  the 
combustion  gas  at  high  temperatures .  Likewise  for  the  oxidizer,  except  blend¬ 
ing  the  oxidizer  vapor  line  with  the  oxidizer-rich  branch. 

EMGML  Input  Data 

Subroutine  ES'GEaL  performs  an  engine  balance  hich  solves  steady-state  pres¬ 
sures  and  flowrates  based  on  best  estimates  of  combustion  efficiency.  It 
also  performs  all  its  own  input  and  output  functions,  except  for  propellant 
properties. 

Types  of  Engine  Balance  Solutions.  Two  types  of  -.ngine  balance  solution^  may 
be  selected  using  indicator  NTYPEB.  Type  1  is  f • c  simulating  an  engine  dtb 
fixed  design  parameters  in  vluch  performance  is  dependent  on  the  pressures  at 
vhich  the  propellants  are  supplied  to  the  feed  systu-m.  Type  2  is  for  a  de¬ 
sign  type  analysis  in  vhich  the  chamber  injector  end  pro  sure  and  mixture 
ratio  of  injected  propellant  flowrates  is  specified,  letting  feed  system  in¬ 
let  pressures  vary  as  the  case  requires. 

Injector  Orifice  Areas.  The  summation  of  injector  orifice  areas  is  required 
by  KS'GHAL  for  both  fuel  and  oxidizer  elements.  Input  specifications  are  set¬ 
up  for  the  program  to  calculate  and  sum  the  areas  of  circular  orifices  from 
diameters  and  uumber  of  orifices.  If  the  iujector  orifices  are  not  circular, 
any  input  combination  which  will  give  the  sane  total  area  is  permissible. 

Operational  Parameters .  Input  of  propellant  valve  inlet  pressures,  PVALVT 
and  PVA1.V0.  is  required  for  type  1  engine  balance  solution  only,  since  these 
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pressures  are  solved  in  EXGRAL  for  a  type  2  solution.  The  mixture  ratio  of 
injected  propellant  flowrates,  XMRI,  and  injector  end  pressure  are  required 
only  for  type  2  solutions,  but  are  used  if  entered  as  estimated  values  for 
type  1  solutions.  The  program  contains  logic  for  predicting  the  values  for 
XMRI  and  PIE;  however,  computational  efficiency  can  be  improved  by  supplying 
better  estimates  of  these  parameters.  The  ratio  of  injector  end-to-nozzle 
stagnation  pressures  is  also  an  optional  input,  which  may  improve  computa¬ 
tional  efficiency  if  a  good  estimate  is  known. 

Efficiency  Factors.  Estimate  values  are  required  for  propellant  vaporization 
efficiencies,  RVAFF  and  RVAP0,  and  for  gas  mixing  efficiency,  ECSMIX,  which 
is  the  ratio  of  mass  flowrate  weighted  theoretical  c*'s  of  individual  stream- 
tubes  with  the  theoretical  c*  at  the  mean  gas  mixture  ratio  at  the  throat. 

If  no  better  estimates  are  available,  use  values  between  0.9  and  1.0.  The 
energy  factor,  ECSENIt,  is  an  input  constant  which  is  to  be  used  to  account 
for  steady-state  losses  such  as  heat  transfer  from  the  combustion  gases  to 
the  chamber  wall.  Also,  this  factor  is  used  for  projecting  losses  when  the 
chamber  vail  is  not  up  to  steady-state  temperature  during  a  duty  cycle  analy¬ 
sis.  If  a  value  is  not  well  known,  ECSENR  might  be  varied  for  correlating 
simulated  performance  with  duty  cycle  test  data. 

Tolerances  on  Solution.  Satisfactory  solutiou  of  an  engine  balance  is  deter¬ 
mined  by  comparing  absolute  differences  in  predicted  and  corrected  values  of 
HVAFF,  HVAPli ,  ECSMIX  and  RPCIN  with  a  specified  tolerance,  EBT0L.  Since  this 
is  not  directly  a  tolerance  on  performance,  it  should  be  smaller  than  the 
accuracy  required  on  the  solution  of  absolute  perf o nuance.  A  tolerance  of 
0.002  vas  used  during  checkout  and  model  evaluation  without  any  difficulty, 
and  this  value  should  yield  results  within  approximately  0.005  (0.5£)  of  a 
fully  converged  solution.  A  second  tolerauce,  EBT0L2,  is  used  as  a  criteria 
for  making  a  second  pass  through  both  the  LISP  model  and  the  single  stream 
tube  analysis  of  PMSTC.  This  may  be  necesr^ry  because  mass  distribution  and 
mean  spray  drop  sizes  are  dependent  on  both  the  fuel  and  oxidizer  injection 
flowrates,  which  ore  in  turn  dependent  on  estimated  vaporization  and  mixing 
efficiencies  and  pressure  ratio.  If  the  corrected  values  of  injection  flow¬ 
rates  deviate  significantly,  more  than  EM0L2,  from  the  predicted  values. 


then  both  the  mass  distribution  and  mean  spray  drop  size  analysis  and  the 
single  stream  tube  combustion  analysis  must  be  updated. 


Fluid  F ] ov  Resistance  Factors.  Steady-state  pressure  losses  due  to  propellant 
fluid  friction  in  the  feed  system  elements  is  calculated  using  the  following 
standard  equation  for  pipe  flow: 
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where  is  the  friction  factor,  which  is  a  required  data  input.  The  program 
user  may  have  a  better  feel  for  the  pressure  loss  Jor  a  reference  flow  con¬ 
dition  than  for  the  friction  factor;  in  which  case,  R^  may  be  back-calculated 
from  the  foreg'  %  equation. 


The  equation  for  the  entrance  pressure  loss  of  the  injector  orifices  is: 
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where  K  in  the  input  entrance  loss  coefficient.  An  input  option  on  K  exists 
for  the  situation  iu  which  the  injector  oriiice  entrance  loss  is  the  primary 
feed  system  pressure  loss  and  the  entire  fuel  (or  oxidizer)  feed  system  pres¬ 
sure  loss, AP  ,  is  known  for  a  reference  flow  conditi  i  A  negative  sign 

i’O 

un  the  input  is  the  indicator  for  the  option,  and  the  absolute  value  of  the 
input  should  be: 
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The  entrance  coefficient  is  then  back-calculated  by  the  program. 
LISP  Input  Data 


Subprogram  LISP  calculates  mass  distribution  and  mean  spray  drop  sires.  Parameters 
for  both  the  injector  design  and  the  mesh  system  for  the  axial  plane  in  which 


the  mass  distribution  is  calculated  (the  "collection"  plane)  must  be  defined 
with  input  data.  Spray  mass  fluxes  emitted  from  injector  elements  are  calcu¬ 
lated  and  summed  to  obtain  the  total  mass  flux  at  each  (r,B)  mesh  point  in 
the  "collection"  plane. 

"Collection"  Plane  Sectors.  In  order  to  minimize  the  number  of  mesh  points 
required  in  the  "collection"  plane  and  the  corresponding  number  of  calcula¬ 
tions,  the  circular  cross-sectional  area  of  the  chamber  at  the  "collection" 
plane  is  split  into  like  sectors  as  small  as  symmetry  conditions  permit. 

Then  calculations  are  made  to  analyze  a  single  sector.  For  plotting,  the 
analysis  is  extended  to  cover  the  complete  circular  area  by  considering  the 
symmetry  as  -  jcified  by  the  JSfM  indicator.  A  mirror  image  type  of  symmetry 
is  specified  if  J3YM=1,  and  a  repeating  image  if  JSYM=2.  For  completing  the 
circle  vith  mirror  image  symmetry,  a  sector  for  repeating  image  symmetry  is 
constructed  first  by  joining  the  analyzed  sector  vith  its  mirror  image.  Then 
the  sector  for  repeating  image  symmetry  is  joined  together  vith  duplicate 
sectors  to  complete  the  circle.  From  these  considerations,  it  follovs  that 
the  defined  sector  must  be  an  even  increment  of  a  half  circle  for  -TSYM=1  and 
an  even  increment  of  a  vhole  circle  for  JSlfM=2. 

Methods  for  Analyzing  the  Sector.  One  of  tvo  methods  must  be  selected  for 
analyzing  a  sector.  The  more  generally  recommended  method  is  to  include  in 
the  analysis  all  injector  elements  vhich  contribute  significantly  to  the  sec¬ 
tor  of  concern,  including  elements  outside  the  angular  range  of  the  sector. 
This  is  the  only  valid  method  for  cases  vhere  elements  are  located  near  the 
axis  of  the  chamber.  The  second  method  is  to  include  only  the  injector  ele¬ 
ments  vithin  the  angular  range  of  the  sector  and  to  "collect"  all  the  mass 
fLvrate  from  those  elements,  even  beyond  the  angular  range  of  the  sector. 

The  mass  flovrate  "collected"  outside  the  sector  is  tn^n  foldev  back  inside 
by  the  program  in  a  manner  compatible  vith  the  type  cf  symmetry  specified. 
Mass  flovrate  collected  outside  a  hard  boundary,  i.e.,  chamber  wall  or  rad¬ 
ial  baffle,  is  folded  into  the  boundary  as  though  the  spray  had  hit  the  vail 
upstream  and  flowed  ulong  the  vail  to  the  collection  plane. 
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Tvo  sectors  are  defined  in  the  input,  one  for  the  "collection"  mesh  system 
and  the  other  for  representing  an  incremental  sector  of  the  chamber  cross  ■ 
section.  These  two  sectors  are  sometimes  referred  to  as  the  outer  and  inner 
sectors,  respectively,  even  though  the  boundaries  may  coincide.  The  "collec¬ 
tion"  sector  mesh  system  contains  NRML  radial  mesh  positions  and  NTUML  ar.gu- 
lar  mesh  positions.  The  first  angular  position  corresponds  with  0=  THETAR. 
Inner  sector  boundaries  are  along  the  WALL  radial  mesh  position  and  along 
the  iSTUR  and  NTHL  angular  mesh  positions.  Generally,  the  outer  sector  ex¬ 
tends  beyond  the  inner  sector  in  the  radial  direction  to  account  for  spray 
hitting  the  wall.  The  radial  boundaries  differ  only  if  the  second  method 
for  analyzing  a  sector  is  selected,  in  which  case  the  flowrate  "collected" 
outside  the  inner  sector  is  folded  back  inside. 

Tht  number  of  mesh  points  in  the  outer  sector  is  the  product  of  NEML  and  NTUML. 
This  product  is  limited  to  400  in  any  combination.  For  CRT  plots,  the  inner 
sector  is  restricted  further:  NWALL<20,  (NTHL-NTHR+l)<20  for  JSYM=1  and 
<39  for  JSYM=2,  and  t!  »•  sector  angle  >5°  for  JSYM=1  and  >10°  for  JSYM-2, 

Spray  Drop  Sizes.  If  the  injector  element  is  specified  as  being  Type  1 
through  Type  5,  LISP  will  calculate  a  mass  median  drop  diameter  for  the  pro¬ 
pellant  of  each  orifice  of  the  element.  These  calculations  are  based  upon 
correlations  derived  from  hot  wax  experiments  with  a  liquid  physical  property 
correction  term,  (fUr/p)  added.  The  coefficient  CDRAR  is  a  multiplier  on 
these  diameters  and  may  be  used  to  orrtlate,  if  necessary,  results  with 
test  data.  Alternatively,  by  assigning  a  value  greater  than  zero  to  the  in¬ 
dicator  IDBAR,  the  program  user  assigns  his  own  estimation  of  drop  diameters. 
For  elements  defined  as  Type  8,  the  user  always  supplies  his  own  estimation 
of  drop  di arte ter.  The  appropriate  mean  droplet  diameter  is  the  mass  median 
diameter. 

Orientation  of  an  Injector  Element's  Coordinate  System.  Individual  (x,y,z) 
coordinate  systems  ure  defined  for  each  element  with  the  origins  located  at 
the  iiapiugenent  points  of  the  elements.  Single  impinging  element  coordinate 
systems  are  specified  in  Fig.  10  .  For  the  doublet  and  triplet  elements,  the 
x  axis  will  correspond  to  the  long  axis  of  the  spray  fan  forced  by  the  im- 
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Figure  10.  Injection  Element  Designations  and 
Coordinate  Systems 


pinging  streams  while  the  y  axis  lies  in  the  plane  of  the  impinging  streams. 
For  the  4-on-l  and  like  doublet  pair  elements,  the  geometries  specified  in 
Fig.  10  were  made  as  consistent  vith  the  doublet  and  triplet  as  was  conven¬ 
iently  possible.  The  definitions  of  x  and  y  axes  for  the  type  8  element  are 
equivalent  to  those  specified  for  the  doublet. 

The  angular  orientation  of  the  (x,y,z)  coordinate  system  of  the  individual 
element  to  the  chamber  coordinate  system  is  defined  in  terms  of  three  rota¬ 
tion  angles  vhose  respective  FORTRAN  names  are  ALFA,  BETA,  and  GAMMA.  If: 

1,  The  element  is  oriented  on  the  injector  such,  that  its  y  axis 
(as  defined  by  Fig. 10  )  is  coincident  vith  the  chamber  radius 
through  the  element  impingement  point,  vith  positive  y  point¬ 
ing  radially  outvard,  and 

2.  The  z  axis  of  the  individual  element  (as  defined  in  Fig. 10  ) 
is  parallel  to  the  thrust  chamber  axis 

then  the  angles  ALFA,  BETA,  and  GAMMA  have  zero  values.  If  the  element  is 
not  oriented  in  this  "basic"  or  reference  alignment,  then 

1.  ALFA  is  the  counterclockvise  angle  of  rotation  around  the 

z  axis  from  its  original  alignment  with  the  chamber  radius. 

2.  BETA  represents  counterclockwise  rotation  around  the  y  axis 
(in  its  uev  positiou  after  the  first  rotation);  finally, 

3.  GAMMA  represents  counterclockwise  rotation  around  the  x 
axis  (in  its  transformed  position  after  the  first  two 
rotations) . 

For  each  rotatiou,  "counterclockv- so"  implies  the  rotation  direction  that 
would  be  seen  by  an  observer  looking  along  the  positive  axis  toward  the  ori¬ 
gin  of  the  clemeut.  Each  element  is  considered  to  consist  of  two  equivalent 
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orifices  designated  as  1  and  2.  The  physical  orifices  which  correspond  to 
these  number  designations  are  shown  in  Fig.,  F10. 

The  LISP  program  also  labels  the  injected  propellant  as  either  1  or  2.  The 
fuel  must  be  chosen  as  propellant  1  and  the  oxidizer  as  propellant  2  because 
of  the  combustion  gas  properties  vs  mixture  ratio  tables  and  because  they  are 
expected  to  bear  these  designations  when  data  are  transferred  from  LISP  to 
PMSTC. 


"Collection"  Plane  Selection.  There  are  no  firm  criteria  for  determining  the 
appropriate  axial  location  for  the  "collection"  plane.  The  distance  down¬ 
stream  of  the  mean  propellant  impingement  point  for  the  collection  plane  is 
generally  selected  within  a  range  of  less  than  l/2  inch  to  2  or  3  inches. 

The  distance  should  be  related  somewhat  proportionately  to  the  physical  size 
of  the  injector  element  and  should  in  general  be  sufficiently  long  enough  to 
allow  the  spray  to  spread  out  across  the  chamber  and/or  to  overlap  with  neigh¬ 
boring  elements.  Also,  the  "collection"  plane  should  be  near  the  region  of 
initial  combustion. 

To  assist  in  evaluating  the  effect  of  "collection"  plane  location,  up  to 
three  locations  may  be  analyzed  in  a  single  computer  run  with  the  use  of 
inputs  ZpM,  Z0M2,  and  Z0M3.  hhen  a  value  of  zero  is  encountered,  the  addi¬ 
tional  LISP  cases  are  not  run.  The  fin?l  case  is  the  one  transferi'ed  to 
the  PMSTC  subprogram  analysis. 

PMSTC  Input  Data 

Input  data  for  the  PMSTC  subprogram  block  are  read  by  subroutine  C INPUT, 

Stream  tube  data  is  initialized  from  LISP-geuerated  mesh  point  flow  data 
via  scratch  data  set  2  if  ILISP/0;  otherwise,  it  is  read  via  punched  cards. 

Size  Control  Integers.  Stream  tubes  formed  from  LISP-generated  data  are 
grouped  into  XOZ0.V  annular  zones  plus  a  wall  boundary  zone  if  NjfZffet  is  less 
than  Ntth’ALL,  The  greater  the  number  of  zones,  the  more  nearly  will  the 


correspondence  of  the  radius  of  an  ensemble  of  spray  in  LISP  be  with  the 
mean  radius  of  the  stream  tube  to  which  it  is  assigned  in  PMSTC.  The  number 
of  zones  does  not  generally  have  a  significant  effect  on  performance.  Each 
N^Z^N  zone  is  divided  into  NSTPZ  stream  tubes  by  grouping  spray  mass  from 
each  LISP  mesh  point  according  to  mixture  ratio.  NSTPZ  should  be  kept  suffi¬ 
ciently  large  to  prevent  too  much  mixture  ratio  averaging,  which  directly 
affects  the  mixing  efficiency.  From  the  standpoint  of  computing  economy,  the 
total  number  of  stream  tubes  should  be  kept  to  a  minimum.  Tbe  program  is 
limited  to  a  maximum  of  19  stream  tubes. 

The  number  of  spray  group  sizes,  NGT,  per  stream  tube  is  limited  to  12,  in¬ 
cluding  both  fuel  and  oxidizer  groups.  If  ILISP^O,  subroutine  SCflMBL  sets 
NGT=12  and  NGF=6,  regardless  of  what  input  values  are  read. 

The  number  of  axial  stations  for  calculations,  starting  at  ZSTART  and  ending 
at  the  throat  plane,  is  specified  with  the  input  variable  NP.  For  read¬ 
ability  of  the  output,  a  whole  number  of  increments  per  inch  is  desirable, 
and  NP  can  be  calculated  simply  as: 

NP  =  ZT-ZSTART 

dZ  1 

vhere  dZ=l/( increment  per  inch)  with  a  typical  range  of  0 ,02<AZt0. 10. 

Combustion  Chamber  Geometry.  The  geometry  of  the  combustor  is  described 
through  the  doubly  subscripted  array  APR0F(J,l).  This  array  is  entered  in 
coordinate  pairs  (L»l,2)  for  each  value  of  J.  The  values  APR0F(J,1)  denote 
axial  distances  from  the  injector  and  the  values  APR0F(J,2)  denote  the 
corresponding  chamber  diameters  at  these  positions.  It  is  required  that 
distance  increase  with  increasing  J  and  that  the  array  progress  from  the 
injector  plane  to  the  nozzle  throat.  That  is,  APR0F(1,1)  is  assumed  to  be 
zero  and  APR0F(NAP,1)  is  the  injector  to  throat  distance.  The  intermediate 
values  with  1>J>NAP  are  used  to  describe  the  wall  profile  of  an  axisymmetric 
chamber.  Cross- sectional  areas  of  z-planes  lying  between  APR0F(J,1)  and  APR0F 
(vM,i)  are  based  on  linear  interpolations  on  diameter.  Therefore,  each  section 
of  the  chamber,  except  the  last,  is  a  section  of  a  cone  or  a  cylinder. 
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The  last  section  pf  the  chamber  is  described  with  a  wall  profile  radius 
of  curvature,  RR,  through  the  throat  region  in  addition  to  the  end 
coordinates,  APR0F(NAP-1,L)  and  APR0F(NAP,L) .  The  resulting  surface,  which 
is  constructed  in  AVAR,  is  a  conical  surface  tangent  with  the  throat  wall 
radius  of  curvature  surface. 

Multiple  stream  tube  PMSTC  analysis  is  continued  past  the  throat  for 

up  to  25  z-increments.  With  the  throat  plane  denoted  by  Z(NP),  mirror-image 

symmetry  is  assumed  such  that  areas  at  Z(NP+1)  and  Z(NP-l)  are  equal,  etc. 

Contraction  Ratio  Tolerance.  A  throat  contraction  ratio  tolerance,  CRT0L, 
is  used  as  a  criteria  for  satisfying  throat  boundary  conditions  in  both  the 
single  and  multiple  stream  tube  analyses.  This  control  becomes  redundant 
when  used  with  the  engine  balance  subroutine,  EN’GBAL,  in  which  pressure  and 
flowrates  are  solved  as  a  function  of  the  throat  area  and  are  compatible  with 
the  throat  area  when  the  predicted  performance  parameters  match  the  corrected 
ones.  For  some  reason,  unknown  at  this  time,  the  contraction  ratio  obtained 
with  the  ENTrBAL  solution  differs  from  unity  by  1  or  2  percent.  Therefore,  an 
adequate  tolerance  must  be  allowed  on  CRT$L  (0.05  recommended)  to  prevent  pro¬ 
gram  termination  due  to  nonconvergence. 

In  the  single  stream  tube  analysis,  the  flow  area  at  the  throat,  computed 
from  gas  continuity  with  the  °ras  velocity  set  equal  to  the  local  sonic  velo¬ 
city,  is  divided  by  the  geometric  throat  area  to  obtain  the  throat  contrac¬ 
tion  ratio.  For  an  acceptable  solution,  the  absolute  deviation  of  this  con¬ 
traction  ratio  from  unity  must  be  less  than  CRT0L ;  otherwise,  parameters  are 
adjusted,  and  a  new  pass  through  the  axial  increments  is  performed. 

For  the  multiple  stream  tube  analysis,  the  throat  flow  area  is  calculated 
differently,  but  the  tolerance  is  still  CRT0L.  The  throat,  or  choked  flow, 
area  here  is  taken  as  the  minimum  cross-secticn  in  the  throat  region.  Also, 
if  the  contraction  ratio  deviates  less  than  three  times  CRT0L,  the  next  pass 
through  the  axial  increments  is  shortened  by  starting  at  the  beginning 
of  convergence. 
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PULSE  Input  Data 


The  pulse  characterization  model,  PULSE,  performs  a  transient  simulation  of 
several  sequences  of  "standard"  width  pulses  and  constructs  performance 
tables  from  which  pulse  performance  of  individual  pulses  can  be  synthesized. 
The  input  data  must  specify  values  for  the  independent  parameters  of  these 
tables,  operating  conditions,  design  data  for  transient  performance  and  em¬ 
pirical  constants. 

Size  Control  Integers  and  Control  Indicators 

The  number  of  pulses  in  each  sequence  of  "standard"  pulse  widths  is  speci¬ 
fied  with  NSPACE.  The  parameter  which  is  varied  between  these  pulses  is  the 
pulse  spacing,  or  electrical  "off-time'  from  one  pulse  to  the  next.  A  maxi¬ 
mum  of  12  pulses  in  a  sequence  is  allowed  fer  covering  a  pulse  spacing  range 
from  the  minimum  spacing  required  up  to  a  spacing  above  which  performance  is 
not  significantly  affected.  Once  the  propellant  in  the  feed  system  has 
boiled  off,  performance  as  a  function  of  pulse  spacing  has  stabilized. 

The  number  of  pulse  sequences  is  specified  with  NTWALL.  Chamber  wall  tem¬ 
perature  is  the  parameter  varied,  hut  only  its  effect  on  heat  soak-back  to 
the  propellant  in  the  feed  system  for  boil-off  is  accounted  for.  Energy  loss 
variations  with  chamber  wall  temperature  is  nut  accounted  for  here,  but  is 

in  the  duty  cycle  analysis.  A  maximum  of  six  pulse  sequences,  or  chamber 
wall  temperatures,  is  allowed,  and  the  temperatures  must  cove-  the  range 
from  ambient  to  steady-stale  wall  temperatures. 

The  modified  Seamaus '  ignition  model,  IGN,  may  be  bypassed  by  entering  a 
”0"  indicator  for  IIGN.  Experience  with  this  model  indicates  that  numerical 
difficulties  will  cause  it  to  fail  and  terminate  prior  to  attaining  ignition 
for  most  cases.  Therefore,  a  "0"  for  XIGN  is  recommended. 

An  indicator  for  the  boil-off  model,  IBdlL,  is  provided.  Without  the  boil- 
off  model,  most  of  the  performance  variation  due  to  pulsing  is  neglected. 
Whereas  this  indicator  was  useful  during  computer  program  checkout,  an 
entry  of  Ml"  for  1D01L  is  recommended  in  order  to  include  the  boil-off 
analysis. 
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Heat  Transfer  Parameters.  The  heat  transfer  considerations  for  pulse-mode 
operation  of  an  attitude  control  engine  is  extremely  complex  and  significant¬ 
ly  affects  pulse  total  impulse  and  mean  specific  impulse  performance.  The 
PMPM  program  does  not  attempt  to  model  the  complex  heat  transfer  conditions; 
but,  instead,  utilizes  gross  heat  transfer  coefficients  and  parameters  which 
must  be  (1)  analyzed  with  a  separate  heat  transfer  analysis  computer  program, 

(2)  be  determined  experimentally  and/or  (3)  be  used  as  correlation  parameters. 
There  are  several  steady-state  heat  transfer  parameters:  a  bulk  chamber  wall 
temperature  (TCWSS),  fuel  and  oxidizer  feed  system  hardware  temperatures 
(TFSBF  §  TFSB0)  and  heat  soak-back  rates  (QSBFSS  for  fuel  and  QSB0SS  for 
oxidizer)  for  heat  transfer  from  the  feed  system  hardware  to  the  propellants 
Chamber  wall  temperature  transients  are  handled  in  the  duty  cycle  analysis 
rather  than  in  the  transient  combustion  analysis.  Initially,  the  duty 
cycle  chamber  wall  temperature  is  at  ambient  temperature,  and  the  wall 
is  assumed  to  heat  up,  while  the  chamber  is  fired,  exponentially  approaching 
steady- stare  wall  temperature.  Similarly,  the  wall  is  assumed  to  cool  while  the 
chamber  is  not  being  fired,  exponentially  approaching  ambient  temperature. 

The  heating  and  cooling  rates  are  controlled  by  the  program  user  through 
the  exponential  coefficients  C0EHTH  and  C0EHTC . 

Ignition  Delay  Constant.  An  ignition  delay  constant,  TAUIGN,  may  be  input 
in  lieu  of  running  the  modified  Seamans'  ignition  model,  IGN.  A  combustible 
gas  mixture  ratio  must  exist  at  the  start  and  during  the  ignition  delay  period 
in  order  to  attain  ignition.  Substantiated  values  of  Ignition  delay  may  be 
difficult  to  obtain.  Typical  values  are  believed  to  be  on  the  order  of  a 
fraction  of  a  millisecond  . 

.Propellant  Valve  Parameters.  Propellant  flow  is  assumed  to  be  controlled  by 
electrically  activated,  solenoid  type  valves.  Valve  opening  and  closing 
response  periods  arc  distinctly  split  between  coil  energize,  or  de-energize, 
time  and  valve  travel  time,  with  the  former  time  consuming  by  far  th*.  major 
portion  of  the  period.  Typically,  for  a  valve  openi ig  response  time  of  5 
milliseconds,  the  valve  travel  time  is  approximately  only  1  millisecond. 

Valve  coil  energize  and  de-energize  response  times  are  calculated  as  linear 
functions  of  applied  voltage,  viU\  the  coefficients  being  part  of  the  input 
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data.  Valve  opening  and  closing  movement  is  specified  as  part  of  the  input 
data  in  the  form  (if  tables  with  fraction  of  restricted  flow  area,  compared 
with  full  open  area,  as  a  function  of  time  from  the  completion  of  coil  ener¬ 
gize  or  de-energize  time.  Considering  the  short  durations,  linear  (two- 
point)  opening  and  closing  travel  functions  are  generally  adequate. 

Propellant  Feed  System  Elements.  The  propellant  feed  system  flow  passages  are 
defined  in  elements:  a  line,  manifold  and  injection  orifices  for  both  the  fuel 
and  oxidizer  branches.  Cross-sectional  area,  length  and  volume  must  be  known 
to  model  steady-state  and  transient  fluid  flow  and  the  propellant  boil-off 
between  pulses.  For  the  line  and  orifice  elements,  diameter  and  length  must  be 
specified  as  input  data,  with  the  volume  implied  assuming  a  constant  cross- 
sectional  area.  For  the  manifold,  however,  the  cross-sectional  area  may  not  be 
constant,  nor  the  flow  path  length  one-dimensional.  Therefore,  volume  and 
cross-sectional  area  at  the  injector  are  the  required  input  data  for  the 
manifold  element  with  an  "effective"  length  implied.  This  length,  relative 
to  the  area,  is  normally  short  enough  to  be  an  insignificant  modeling 
parameter. 

DCYCLE  Input  Data 

« 

A  mission,  or  duty  cycle,  for  a  pulse-mode  engine  must  be  fully  specified 
with  input  data  describing  sequentially  pulse  widths  (electrical  on-times) 
and  spacings  (electrical  off-times). 

Pulse  Sequences.  Frequently,  a  pulse-mode  engine  is  fired  in  bursts  of  con¬ 
stant  width,  equally  spaced  pulses.  A  "burst"  of  pulses  of  this  type  is  re¬ 
ferred  to  here  as  a  "sequence"  of  pulses.  To  minimize  input  data,  duty  cycles 
are  specified  in  terms  of  sequences,  where  the  number  of  pulses  in  a  sequence 
may  be  as  few  as  one.  The  input  variable,  NSEQ,  specifies  the  number  of  pulse 
sequences  iw  the  duty  cycle.  Each  sequence  is  described  with  the  number  of 
pulses  (NTS),  pulse  width  (iVIDTEl),  pulse  off-time  within  the  sequence  (j)?FB) , 
and  pulse  off-time  between  the  la3t  pulse  and  the  first  pulse  of  the  next  se¬ 
quence  (^FFC).  For  the  last  sequence,  $FFC  should  be  set  equal  to  the  maxi¬ 
mum  off-time  in  the  PSPACE  array  in  subprogram  block  PULSE. 
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PROGRAM  OPERATING  INSTRUCTIONS 


Program  operating  instructions  are  presented  in  this  section  to  specify 
specific  conditions,  other  than  input  data  instructions,  required  in  pro¬ 
cessing  a  computer  run  with  the  PMPM  computer  program. 

DECK  SETUP 

Because  of  the  very  large  size  of  the  PMK-1  computer  program,  the  program  is 
often  executed  in  an  overlay  mode  or  in  separate  parts.  This  may  be  done 
for  several  reasons:  (l)  to  stay  vithin  the  computer  core  size  limitation, 
(2)  to  reduce  core  size  usage  with  overlay  when  the  computer  is  run  in  a 
MVT  (multiple  variable  task)  environment,  or  (3)  to  reduce  deck  size  vith  a 
partial  program  for  ease  of  handling. 

Overlay  Structure 


PMPM  overlay  structure  is  shown  in  Fig.  11  .  The  horizontal  lines  show  the 
level  of  the  overlay  vith  each  branch  set  up  to  occupy  core  separately  in 
time  from  the  other  branches  attached  to  the  same  horizontal  line.  Each 
subprogram  or  label  common  block  name  (which  are  enclosed  in  slashes)  occurs 
only  in  one  segment.  Segments  are  numbered  in  circles  on  the  overlay  chart. 

The  root  segment,  segment  1,  contains,  in  general,  the  main  control  program 
and  utility  programs  which  are  used  extensively  in  many  brunches  of  the  over¬ 
lay  structure.  The  primary  branches,  which  connect  at  the  overlay  A  level, are 
the  LISP  model  in  segments  2  through  5,  the  PMSTC  model  in  segments  6  through 
11,  and  the  PULSE  and  DCYCLE  models  in  segments  14  through  20. 

The  recommended  method  for  specifying  the  overlay  structure  when  executing 
the  program  (with  TDK  excluded)  is  to  include  in  the  dock  setup  the  small 
deck  of  overlay'  and  insert  curds  (which  are  provided)  at  the  end  of  the 
link  edit  section.  Insert  cards  are  not  provided  for  the  TDK  program  block, 
which  contains  approximately  100  subroutines  and  many  label  common  blocks. 
Instead,  the  order  of  the  delivered  deck  is  arranged  to  correspond  with  the 
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Figure  11  .  PMPM  Program  Overlay  Structure 
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Figure  11.  PMPM  Program  Overlay  Structure  (Continued) 


TDK  overlay  structure.  The  overlay  cards  for  running  TDK  separately  from 
FMPM  are  to  he  inserted  in  the  object  deck  at  the  head  of  each  overlay 
branch;  i.e.,  with  2  overlay  A  cards,  5  overlay  B  cards,  and  2  each  overlay 
C  and  D  cards.  ■ 

The  overlay  structure  was  setup  to  minimize  the  computer  core  size  require¬ 
ments  without  incurring  excessive  overlaying  during  execution.  Except  for 
the  TDK  program  block,  the  segments  of  the  program  are  generally  loaded  se¬ 
quentially  with  a  minimum  of  subsequent  reloading.  Overlay  is  not  used  in 
iterative  or  incremental  step- type  solutions.  However,  overlay  efficiency 
on  computers  having  more  than  the  minimum  required  core  capacity  can  be  im¬ 
proved  by  consolidating  the  segments  to  just  fill  the  available  core  size. 

Partial  Program  Configurations 

The  generally  recommended  "complete"  PMPM  computer  program  configuration  ex¬ 
cludes  (vith  the  use  of  dummy  subprograms)  the  very  large  TDK  subprogram 
block  in  PMSTC  and  the  IGN  model  in  PULSE,  since  the  nozzle  performance  anal¬ 
ysis  vith  TDK  is  normally  either  executed  separately  or  omitted  and  since  the 

* 

ignition  analysis  with  IGN  is  inoperable  for  the  general  case. 

i’he  computer  core  size  required  for  this  "complete"  configuration  without 
overlay  is  approximately  115K  words,  where  the  unit  K  is  1024.  The  required 
core  sizes  quoted  here  can  only  oe  approximate,  since  size  is  variable  in 
respect  to  compiler  efficiency,  computer  system  software  and  buffer  sizes, 
and  the  word  structure.  Sizes  quoted  here  are  taken  from  operational  data 
run  on  an  IBM  System  360,  model  GS  computer  system  vith  the  units  of  bytes 
converted  to  words  by  dividing  by  four.  The  same  configuration  vith  overlay 
requires  approximately  53K  words  of  core. 

A  natural  divis’ m  of  the  FMPM  computer  program  is  between  the  steady-state 
model,  PMDE31,  and  the  pulse  characterization,  PULSE,  and  duty  cycle,  DCYCLE, 
models.  As  a  rule,  the  steady-state  results  should  be  reviewed  before  pro¬ 
ceeding  into  the  pulse  analysis,  in  which  case  the  program  deck  might  just 
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as  veil  be  dividecj.  to  minimize  the  deck  size,  and  the  number  of  overlay  seg¬ 
ments  and/or  required  core  size.  With  dummy  PULSE  and  DCYCLE  subprogram 
blocks,  the  core  requirement  without  overlay  is  approximately  90K  words. 

With  a  dummy  PMDHt  subprogram  block,  the  core  requirement  without  overlay 
is  approximately  66K  vords. 

The  CRT  plotting  routines  comprise  a  sizable  section  of  the  source  program. 
These  routines  are:  LINEG,  LGRID,  and  CRTIiBL  in  segment  1  of  the  overlay 
structure,  all  of  them  in  segment  5,  STCRT  in  segment  8,  PT (/TT  in  segment 
10,  and  PULPLT  in  segment  19.  Also,  there  is  an  extensive  number  cf  library 
routines  referenced  by  these  supplied  routines.  If  CRT  plotting  is  net 
supported  by  the  user's  computer  facility,  the  CRT  routines  should  be  removed 
from  the  program  deck  and  replaced  with  the  following  subroutines,  either 
dummy  or  rewritten  to  be  compatible  with  the  user's  plotting  equipment: 

PL^TC ,  PL/TN,  PL/T2,  STCRT,  PL/TT,  and  PULPLT.  Also,  the  plotting  routines 
referenced  on  the  overlay  insert  cards  must  be  deleted  or  supplemented  as 
the  modifications  require. 

DATA  SET  USAGE 

» 

For  normal  I/O  (input/output)  operatious  in  the  FMPM  computer  program,  the 
standard  data  set  reference  numbers  5  and  6  are  used  with  the  transmittal 
statements  READ  and  WRITE,  respectively.  These  data  sets  may  refer  to  se¬ 
parate  magnetic  tape  drive  mits  or  a  field  on  a  disk  pack.  The  transmittal 
statement  PUNCH,  for  which  the  computer  system  designates  the  data  set  num¬ 
ber,  is  used  for  the  punched  card  output.  (At  Rocketdyne,  data  set  14  is 
used  for  punched  card  output,  whereas  data  set  7  is  more  generally  used  for 
this  function.)  If  the  CRT  plotting  routines  are  used,  another  data  set 
(number  16  at  P.ouketdyne)  is  required,  which  is  also  system  designated. 

Additional  data  sets  are  required  to  transfer  data  between  subprograms  and 
to  store  and  accumulate  large  quantities  of  data  fox-  delayed  processing. 

Those  special  data  sets,  which  must  be  defined  when  running  PJIPM,  are  x-e- 
ferenccd  as  numbers  2,  3,  and  4  when  runuing  without  TDK.  With  TDK,  10  and 
11  ore  also  required.  Table  3  describes  special  data  set  usage,  the  sub¬ 
programs  where  the  data  sets  are  referenced  and  variable  names  used  in  re¬ 
ferencing  data  sets. 
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TABLE  3  .  SPECIAL  DATA  SET  USAGE 


Data  Set 

No. 

Var.  Name 

Usage 

2 

M 

(1)  Transfer  data  from  LISP  to  STAPE  in 

PMSTC. 

(2)  Save  data  in  PMSTC  for  plotting  in  STCRT , 

KSTC 

(3)  Transfer  data  from  TDKSET  in  PMSTC  to 

TDK  subroutines  BRAIN,  LINE,  and  CHAR. 

3 

- 

(l)  Subroutine  ITER8  saves  and  retrieves 
data  for  iterating  in  PMSTC. 

KPSS 

(2)  Transfer  data  for  PSS  array  in  TDK  from 
subroutine  BRAIN  to  LINE. 

4 

JTAPE 

(1)  It.  PMSTC,  CPM1  saves  droplet  mass  flow¬ 
rates  and  residence  times  and,  in  PMDER, 
CPM2  uses  this  data  for  calculating  the 
spray  depletion  functions. 

(2)  Subroutine  TTAPE  in  TDK  stores  master 
thermodynamic  data  from  block  data  and 
retrieves  it  in  subroutine  STFS'ET. 

10 

KREAX 

Transfer  reaction  rate  tables  in  TDK  from 
subroutine  INPUT  to  PACK1P. 

11 

KSTF 

Transfer  specific  thermodynamic  data  from 

JANAF  unit  in  TDK  from  STFSLT  to  TUERM0, 

i _ 

PACK1P  and  THERM. 

PROGRAM  EXECUTION  LIMITS 


Limits  are  generally  specified  on  several  computer  operating  functions  to 
terminate  computer  execution  for  nonstandard  or  excessive  operation.  In 
particular,  if  the  program  execution  gets  in  an  endless  loop  for  some  reason, 
the  limits  will  terminate  execution  instead  of  letting  the  computer  get  hung¬ 
up.  Limits  which  are  normally  specified  include  execution  time,  printed  out¬ 
put  line  count,  and  number  of  CRT  frames.  Normally  limits  are  specified  with 
considerable  margin  to  prevent  termination  caused  by  a  miscalculation. 

Limits  vary  so  greatly  with  specific  cases  run  that  only  rough  estimates  can 
generally  be  made  from  sample  cases.  For  the  sample  case  shown  in  Volume  III 
of  the  program  documentation,  the  number  of  lines  of  printed  output  is  approx¬ 
imately  18,000  and  the  number  of  CRT  frames  was  30.  The  CRT  output  can  be 
determined  precisely  from  the  input  data,  but  a  margin  for  miscalculation  is 
recommended  here,  also. 

The  limit  on  execution  time  at  Roclcetdyne  is  based  on  CPU  time,  not  clock 
time,  since  the  computer  is  run  in  an  MVT  environment.  Computer  usage  for 
cost  allocation  at  Rocketdyne  is  measured  in  "billing  units"  (BU's).  These 
BU's  are  calculated  with  a  formula  which  includes  CPU  and  channel  times, 
core  region  size  and  amount  of  peripheral  equipment  used.  A  BU  is  roughly 
equivalent  to  a  minute  of  clock  time  for  a  standard  engineering  job  run  by 
itself  on  an  IBM  system  360,  model  65  computer.  The  documented  sample  case 
used  11.5  BU's,  and  the  CPU  time  was  approximately  8  minutes.  A  typical 
TDK  case  requires  40  or  more  BU's,  and  an  IGN  case  requires  on  the  order  of 
2  or  3  BU's. 
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,:  PROGRAM  OUTPUT 

The  PMPM  computer  output  consists  of  printout,  punched  cards  and  CRT 
graphi  al  plots.  Appendix  D,  Volume  III,  contains  printout  from  an  example 
case  which  was  run  using  the  input  data  listed  in  Appendix  C,  Volume  II. 
Card  images  of  the  punched  output  appears  in  the  printout.  Most  of  the 
printout  is  self-explanatory,  but  the  parameters  which  are  labeled  only 
with  their  FORTRAN  code  name  need  to  be  described.  Input  data  are  printed 
out  in  this  manner,  using  the  coded  names  defined  in  Table  1.  Coded  para¬ 
meter  names  in  the  printout  which  are  not  part  of  the  input  data  are 
described  in  Table  4.  A  brief  discussion  of  the  output  by  major  subpro¬ 
gram  blocks  is  presented  in  this  section. 


COPE 

A1IL 

A1SL 

CD 

CFSF.CFS0 

C0MB 

CSTAR 

ECSTAR 
ECSMXE 
EPS1 ,  2 , . . ,  6 
EM 

ERRMR 

ESPIMP 

F 

FF,F0 

IP 

ISEQ 
ITER 
KPRIME 
LI,  L2 

0FF1.0FF2 

0XFRP 
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..... 
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TABLE  4. 

FORTRAN  CODES  OF  PARAMETERS  IN'  PRINTOUT 
WHICH  ARE  NOT  IN  INPUT  DATA  TABLE 

DESCRIPTION 

Slope  of  TDK  initial  line  at  intersection  with  streamline 
Slope  of  streamline  at  intersection  with  TDK  initial  line 
Nozzle  discharge  coefficient 

Fuel  and  oxidizer  feed  system  flowrate  coefficients 
in  ENGBAL 

Logical  variable  for  combustion  (true  or  false) 

Characteristic  velocity,  c*  (used  as  heading  for  c*  table 
and  for  final  %'alue) 

Characteristic  velocity  efficiency  (^c*) 

Estimated  mixing  limited  Dc* 

Area  expansion  area  ratios  (e^,  e^,  .  ..e^.) 

Rupe  mixing  efficiency  factor 

Deviation  of  ENGBAL  mixture  ratio^  predicted  minus 
corrected  values 

Pulse  mean  specific  impulse  efficiency, 

False  s 

Fuel  and  oxidizer  mass  flowrate  continuity  correction 
factors  going  from  LISP  to  PMSTC 

Pulse  number  in  sequence 

F iJ  ??  sequence  number  in  duty  cycle 

Iteration  count  in  ENGBAL 

Table  values  of  static  droplet  evaporation  coefficients 

TDK  initial  line  coordinates  (sequenced  from  nozzle  wall) 
which  bracket  axial,  2,  location 

Electrical  off-time  before  and  after  a  pulse 

•  ■* 

Reduced  oxidizer  fraction,  w/(w  ♦  c  w„),  where  c  is 

o^o  m  m 

the  midpoint  mixture  ratio  in  the  combustion  gas  tables 


TABLE  4  (Coutinued) 


CODE 


DESCRIPTION 


PIEE 

PMF.PM0 

PNS 

P0VER 

PWIDTH 

R1,R2 

RDSL,RDSL2 

RH0F,RH00 

RVAPFE, 

RVAP0E 

SIMEAN 

SPIMP 

SUMSPF, 

SUMSP0 

SUMTI 

SUMW1 , SUMW2 
SUMWF.SUMW0 
T 

TCSTR 

TDR0P 

TGAS 

TI 

TIA.TIB 


Estimated  injector  end  chamber  pressure 

Fuel  and  oxidizer  feed  system  manifold  pressures 

Stagnation  pressure  at  chamber  nozzle 

Maximum  pulse  pressure  overshoot  in  percent  of 
steady- state  pressure 

Pulse  electrical  on-time 

Factors  used  to  corrected  fuel  and  oxidizer  injection  flow' 
rates  in  LISP  to  correspond  with  ENGBAL  flowrates 

Radii  of  a  streamline  (dividing  two  stream  tubes) 

at  ends  of  Z  increment  bracketing  the  TDK  initial  line 

Fuel  and  oxidizer  densities 

Estimated  fuel  and  oxidizer  vaporization  efficiencies 

Cumulative  mean  specific  impulse  of  pulses  in  a  duty  cycle 
Pulse  mean  specific  imputse 

Instantaneous  total  fuel  and  total  oxidizer  mass  weights 
of  spray  in  combustion  chamber 

Cumulative  total  impulse  of  pulses  in  a  duty  cycle 

Total  fuel  and  total  oxidizer  spray  mass  flowrates  in  LISP 

Cumulative  fuel  flow  and  oxidizer  flow  during  a  duty  cycle 

True 

Theoretical  c* 

Pulse  response  time  for  chamber  pressure  to  drop  to  10 
percent  of  steady-state  from  the  off-signal 

Instantaneous  chamber  gas  temperature 

Pulse  total  impulse 

Pulse  start  total  impulse  and  decay  total  impulse 
interpolated  from  tables  characterizing  a  standard 
width  pulse 
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TABLE  4  (Continued) 


CODE 

TIDR0? 

TIF1,TIF2 

TIFACT 

TIRISE 

TMFLF,TMFL0 

T0XFRP 

TRISE 

TRISE2 

TSTART 

TWALL1,2,3 

TWALLM 

VENSF.VENS0 

VFSF,VFS0 

VIL 

WENSF , WENS0 

KF,K0 

KGEXH 


DESCRIPTION 

Total  impulse  during  a  pulse  pressure  drop  response 
period,  TDR0P 

Total  fuel  and  total  oxidizer  injected  into 
chamber  in  LISP 

Pulse  total  impulse  adjustment  factor  to  account  for 
variation  in  combustion  gas  energy  losses  during 
a  duty  cycle 

Total  impulse  during  a  pulse  pressure  rise  period,  TRISE 

Total  fuel  and  total  oxidizer  steady-state  flowrates 

Table  of  reduced  oxidizer  fraction  Tsee  0XFRP)  corresponding 
with  array  of  mixture  ratio,  TMR,  in  input  tables  of 
combustion  gas  properties 

Pulse  response  time  for  chamber  pressure  to  rise  to 
90  percent  of  steady-state  from  the  on-signal 

Pulse  response  time  for  chamber  pressure  to  rise 
from  10  to  90  percent  of  steady-state 

Time  into  a  duty  cycle  of  the  on-signal  for  a  pulse 

Chamber  wall  temperature  during  a  duty  cycle  at  the  start 
of  a  pulse,  at  the  cut-off  of  a  pulse  and  at  the 
end  of  the  pulse  tail-off  period 

Chamber  wall  temperature  at  the  time  half  way 
through  the  on-time  period. 

Velocities  of  fuel  and  oxidizer  spray  ensembles  formed 
during  an  incremental  time  period  in  PULSE 

Total  volumes  in  fuel  and  oxidizer  feed  systems 

Stream  tube  gas  velocity  at  the  S  location  where  a 
dividing  streamline  intersects  the  TDK  initial  line 

Mass  weights  of  fuel  and  oxidizer  spray  ensembles 
formed  during  an  incremental  time  period  in  PULSE 

Total  fuel  and  oxidizer  mass  weight  flows  during  a  pulse 
in  DCYCLE 

Mass  weight  of  gas  flowing  through  the  chamber  nozzle 
during  an  incremental  time  period  during  pulse 
transient  analysis 


TABLE  4  (Continued) 


CODS 

WGF,WG0 

KGFCUM, 

NG0CUM 

X1D 

XMRIL 

XMRMN 

Z,ZL2 


DESCRIPTION 

\  . . . 

Fuel  and  oxidizer  gas  flowrates  from  LISP  to  PMSTC 

Instantaneous  residuals  of  fuel  and  oxidizer  gases  in 
combustion  chamber  during  pulse  transient  analysis 

Axial  chamber  location  for  beginning  of  the  PMSTC  analysis 

Stream  tube  gas  mixture  ratio  at  the  location  where  a 
dividing  streamline  intersects  the  TDK  initial  line 

Cumulative  mean  mixture  ratio  of  propellants  used 
during  a  pulse  duty  cycle 

As  used  in  calculating  TDK  initial  line  parameters, 

1  is  the  first  axial  step  downstream  of  the  intersection 
of  a  streamline  with  the  initial  line  and  ZL2  is  the 
axial  coordinate  of  the  initial  line  at  Z 


ENGBAL  Output 


The  engine  balance  routine,  ENGBAL,  is  generally  called  several  times 
throughout  a  computer  run.  On  the  initial  cali,  the  ENGBAL  input  data, 
coefficients  and  initial  values  are  printed  out.  Injected  mixture  ratio 
is  iterated  on  in  solving  the  engine  balance,  and  a  line  of  printout  is 
generated  for  each  iteration  step.  Each  engine  balance  solution  is 
printed  out.  The  solution  is  based  on  estimated  values  of  spray  vapor¬ 
ization  efficiencies,  mixing  limited  c*  efficiency  and  the  ratio  of  in¬ 
jector  end-to-nozzle  stagnation  pressures.  At  the  end  of  each  PMSTC  run 
estimated  and  calculated  values  of  these  parameters  are  printed  out. 

LISP  Output 

The  output  of  the  LISP  computer  program  is  provided  in  beth  the  form  of 
tabular  printout  and  of  computer-plotted  CRT  graphs.  Also,  data  are 
written  on  scratch  data  sets  to  be  read  and  used  by  the  STC  computer 
program. 

LISP  tabular  output  is  presented  in  the  example  case  printou*.  First, 
there  is  a  tabulation  of  all  input  data  which  permits  both  a  full  docu¬ 
mentation  of  the  computer  run  conditions  for  later  analysis  and  a  convenient 
method  to  check  input  for  errors  if  unusual  results  are  calculated.  The 
input  data  table  is  follrwed  by  a  one  page  table  of  specific  spray  dis¬ 
tribution  coefficient  val"is  used  for  e-ch  element  specification.  As  in 
the  example  ca  e,  this  table  of  coefficients  may  be  preceded  by  a  warning 
message  that  some  unlike  doublet  design  or  flow  parameter  falls  outside 
the  correlation  range  in  the  LISP  library  of  coefficients  along  with 
values  of  pertinent  parameters.  None  of  the  other  coefficient  subroutines 
provides  such  a  message.  Next,  c  table  cross  references  (by  injector 
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element,'  the  calculated  flow  rates  and  drop  sizes  before  evaporation 
w* mv  the  read-in  element  coordinates. 

Following  the  element  reference  table,  there  are  two  extensive  tables 
referenced  to  the  combustion  zone  mesh  points.  The  first  of  these  tables 
lists  the  coordinates  of  the  mesh  points  in  the  chamber  slice  at  the 
collection  plane  zQ,  together  with  the  weight  flux,  the  total  collected 
mass,  the  three  mean  droplet  velocity  components,  and  the  mean  drop 
diameters  of  each  propellant  at  the  mesh  points.  The  total  collected 
mass  at  a  mesh  point  is  defined  as  the  weight  flux  times  the  associated 
area  at  the  mesh  point.  The  values  in  the  first  of  these  tables  are  based 
upon  cold  flow  conditions,  i.e.,  no  vaporization  is  assumed  between 
elements  and  mesh  point.  The  mesh  points  are  listed  in  ascending  order 
according  to  radial  and  angular  coordinate.  The  last  column  of  this  table 
lists  the  sum  of  the  collected  mass  of  each  propellant  at  all  the  mesh- 
points  of  constant  radial  coordinate,  i.e.,  lists  the  radial  distribution 
of  the  spray  mass  flux.  The  second  of  these  tables  again  lists  the 
coordinates  of  the  mesh  points  in  the  chamber  slice  in  the  plane  zQ,  to¬ 
gether  with  the  reduced  weight  fluxes  and  droplet  diameters  of  the  collected 
spray  after  evaporation.  A  mass-weighted  average  evaporation  of  the  original 
spray  flux  to  each  mesh  point  is  also  listed  in  this  table.  At  the  bottom 
of  this  table  are  listed  the  Rupe  mixing  factor,  E  ,  the  mixing  limited 
c*  efficiency  and  the  overall  pex-cent  vaporization  of  each  propellant. 

Samples  of  the  LISP  graphical  output  are  shown  in  Fig.  12  through  16. 

Figure  12  shows  the  mesh  system  for  the  chamber  slice  analyzed  and  the 
element  origin  locations  £o-  all  injection  elements*  considered  to  contribute 
flux  to  that  slice.  Figure  15  is  an  example  of  the  fuel  and  oxidizer  mass 


*  This  case  has  only  a  single  injection  element  which  is  at  the  center. 
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F'OMrc  13.  fuel  and  Oxidiser  Mass  Flux  Profiles  Compjted 
by  LISP  at  a  Given  Chamber  Radius 
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profiles  around  tile  chamber  slice  at  one  fixed  chamber  radius.  Figures  14 
and  15  are  the  contour  plots  of  fuel  and  oxidizer  mass  flux  for  the  entire 
chamber  cross  section.  A  similar  plot  for  total  mass  flux  is  not  shown. 
Figure  16  is  a  contour  plot  of  a  modified  fuel  fraction  function.  The 
expression  plotted  is  given  at  the  top  of  the  figure;  it  was  chosen  be¬ 
cause  it  is  bounded  between  zero  and  unity  and  has  a  values  of  0.5  at 
the  injection  mixture  ratio. 

PMSTC  Output 

A  sample  case  of  PMSTC  computer  program  printout  is  inc  uded  in  Volume  III. 
Input  data  are  written  out  immediately  as  they  are  read  in.  This  documents 
the  data  used  for  the  particular  case  as  well  as  showing  whether  or  not 
the  data  were  read- zn  properly.  The  input  section  should  be  examined 
for  each  case  run  to  be  sure  that  the  intended  input  data  were  actually 
used. 


Input  data  transferred  from  LISP  are  not  printed  out,  but  a  table  of 
diagnostic  data  from  subroutines  STAPE  and  SCRMBL  is  printed.  Parameters 
appearing  in  that  table  are: 

SUMW1,2  -  Total  fuel,  oxidizer  spray  flowrates  summed  over  all 
mesh  points  in  STAPE 

WGF,0  -  Gaseous  fuel,  oxidizer  flowrates 

TIF1,2  -  Total  fuel,  oxidizer  flowrates  transferred  from  LISP 

FF,P0  -  The  ratio  TIF1/(SUMW1  +  NGF) ,  etc. 

K  -  The  number  of  circular  rings  of  mesh  points  from  LISP 

SMBL  -  Total  flowrate  assigned  to  the  wall  boundary  layer 
stream  tube 

PCTT  -  The  ratio  SM6L/(TIF1  ♦  TIF2) 

0TT  -  The  product  PCTBL*(TIF1  ♦  TIF2) 

N1,N2  -  The  indices  of  circular  rings  of  mesh  points  in  a 

given  geometric  zone 

SIW  -  The  cumulative  total  flowrate  in  a  geometric  ion©  and 
those  set  up  prior  to  it 

The  stream  tube  initialization  data  are  tabulated  and  simultaneously  punched 
out  in  cards. 
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Based  on  the  stream  tube  initialization  data,  a  table  is  printed  out  of 
stream  tube  total  flowrates  and  overall  mixture  ratios.  This  table  is 
followed  by  values  of  the  Rupe  mixing  efficiency  factor,  E  ,  and  a 
mixing  c*  efficiency  for  the  stream  tube  flows.  The  latter  represents 
an  upper  limit  for  multiple  stream  tube  c*  efficiency,  since  it 
corresponds  to  complete  evaporation  and  burning  of  sprays  within  all  stream 
tubes.  These  mixing  efficiency  factors  will  differ  from  those  calculated 
in  LISP  because  of  the  gas  distribution  in  PMSTC  and  the  combining  of 
stream  tubes  versus  mesh  point  flows. 

Subroutine  AVAR  sets  up  the  array  of  chamber  areas  and  writes  out  a  table 
of  chamber  geometry  information.  Similarly,  subroutine  KPRIME  computes 
and  writes  out  tables  of  evaporation  coefficients. 

Single  stream  tube  analysis  is  preceded  by  writing  out  a  one-page  table 
of  input  total  flows  and  averaged  spray  and  gas  parameters.  During  single 
stream  tube  analysis,  data  are  written  out  as  they  are  generated.  At 
each  z-plane  to  be  printed,  complete  gas  and  propellant  spray  group  data 
are  giver..  Additionally,  the  percentages  of  propellants  evaporated  and 
burned  are  listed  and  volume  number  mean  propellant  droplet  diameters, 

Djq,  are  computed. 

Two  values  each  of  flow  area  and  contraction  ratio  may  be  given.  Where 
the  gas  flow  is  subsonic,  these  should  agree  with  each  other  precisely, 
whereupon  the  second  set  is  not  printed.  At  or  near  the  nozzle  throat 
plane,  the  two  sets  may  disagree  because  the  gas  velocity  has  been  set 
equal  to  sound  speed  and  the  local  nozzle  area  adjusted  to  satisfy  mass 
flow  continuity.  The  contraction  ratio  calculated  from  continuity  at 
the  throat  plane  is  used  as  a  multiplier  (if  it  differs  from  unity  by 
more  than  CRT0L)  to  adjust  initial  plane  chamber  pressure  for  a  next 
iteration  of  single  stream  tube  analysis.  After  the  throat  plane  data 
are  written  out,  an  engine  balance  is  calculated  and  printed. 
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li/hen  the  foregoing  alnalysis  has  converged  on  its  solution,  a  performance 
summary  sheet  with  c*,  thrust  and  specific  impulse  is  printed.  Next,  the 
input  value  of  nozzle  radius  ratio  and  calculated  value  of  mean  nozzle 
expansion  coefficient,  \  ,  are  used  by  TRANS  to  generate  transonic  flow 
region  isobars.  The  reduced  corrdinates  and.  flow  directions  for  each 
of  20  points  along  each  isobar  are  written  out,  beginning  with  the 
furthest  downstream  isobar  and  progressing  upstream.  Additionally,  for 
thea=0  isobar,  the  absolute  coordinates  are  written  out  for  40  points. 

TRANS  also  generates  a  CRT  plot  of  isobars,  which  is  displayed  in  Fig.  17. 

Multiple  stream  tube  analysis  follows  the  foregoing  single  stream  tube 
and  TRANS  analyses.  Stream  tube  input  data  are  re- initialized  and  some 
additional  data  are  written  out  to  more  completely  define  the  initial- 
plane  conditions.  Initial-plane  pressure  is  taken  from  the  engine  balance 
solution  which  is  performed  just  prior  to  the  multiple  stream  tube  analysis. 

At  each  prescribed  z-plane  for  printing  multiple  stream  tube  results, 
complete  definitive  data  for  combustion  gases  and  propellant  sprays  are 
written  out.  Local  chamber  area  ar.d  contraction  ratio  are  given;  additionally, 
overall  precentages  of  the  propellants  evaporated  and  burned  are  listed. 

At  the  throat  position  and  intermittently  downstream,  diagnostic -type 
printouts  containing  data  concerning  dividing  streamline  intersections 
with  the  a  =0  isobar  are  inserted  between  the  regular  z-plane  printouts. 

A  summary  table  of  these  data  is  given  near  the  end  of  the  multiple  stream 
tube  printout.  Finally,  a  long  summary  table  is  given  of  the  stream  tubes' 
outer  radii  at  each  z-plano.  This  is  terminated  with  the  minimum  value 
of  the  sum  of  stream  tube  areas,  the  ratio  of  that  value  to  throat  area 
ratio  and  the  instantaneous  weight  of  the  gas  in  the  combustion  chamber. 

Anoths.  engine  balance  is  performed  using  the  latest  calculated  values 
of  vaporization  and  mixing  efficiencies,  and  ratio  of  injector  end- 
to-nozsle  stagnation  pressure.  The  magnitudes  of  the  deviations  in  these 
calculated  and  the  previous  predicted  values  determine  whether  or  not 
all,  or  a  portion,  of  STC's  multiple  stream  tube  analysis  will  be  repeated. 

If  so.  it  is  readily  apparent  in  the  printout.  Figure  IS  is  a  computer 
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Figure  17.  Isobars  Showing  No::le  Pressure  Distributions 
Calculate-  c>  Subprograms  TRANS 
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AXIAL  DISTANCE 


Figure  18.  Computer-Plotted  Stream  Tubes  Calculated  in  PMSTC 


plotted  graph  showing  the  stream  tubes’  outer  radii  along  the  entire 
chamber  length. 

Following  the  last  pass  through  the  multiple  stream  tube  analysis,  an 
updated  performance  summary  sheet  is  printed  and  a  table  is  printed  which 
lists  the  data  punched  by  STC  for  subsequent  use  in  running  the  TDK  sub¬ 
program. 

Fuel  and  oxidizer  spray  depletion  functions,  which  are  the  mean  time 
traces  of  the  mass  of  a  spray  droplet,  are  calculated  in  PMSTC  for  use 
in  PULSE.  This  data  are  tabulated  in  the  printout,  punched  on  cards  and 
plotted  by  the  computer  (F;g.  19). 

PULSE  Output 

Following  the  printout  of  the  PULSE  input  data,  an  extensive  table  of 
transient  pulse  performance  data  is  printed  out  (optional)  with  a  one 
page  summary  of  pulse  performance  at  the  end  of  each  pulse.  In  the 
transient  table,  the  following  parameters  are  tabulated:  time,  thrust, 
chamber  pressure,  reduced  oxidizer  fraction  of  chamber  gas,  weights  and 
velocities  of  fuel  and  oxidizer  ensembles  formed  during  time  increment, 
nozzle  flow  for  time  increment,  accumulated  chamber  gases  and  propellant 
•spray,  and  gas  temperature.  Logical  parameter  C0MB  is  also  tabulated, 
which  indicates  whether  combustion  is  present  or  not.  Additional  diag¬ 
nostic  type  data  may  be  printed  out  by  option.  Computer  plots  are  shown 
of  reduced  pulse  thrust  and  pressure  (Fig.  20)  and  flowrates  and  oxidizer 
fraction  (Fig.  21). 

Following  the  printout  of  the  last  pulse,  a  parametric  tab!?  of  pulse 
performance  characteristics  is  constructed  from  the  transient  data.  A 
summary  of  the  table  is  printed  out,  and  the  complete  table  is  punched 
on  cards,  whch  is  also  listed  in  the  printout. 
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Figure  20.  Reduced  Pulse  Thrust  and  Pressure  Calculated 
in  Subprogram  PULSC 
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REDUCED  INJECTION  FLO  IATES  AND  OX  ID  FRACTION 
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Fig  ire  21.  Reduce  Flowrates  and  Oxidiacr  Fraction  Calculated 
in  Suberograa  PULSE 


DCVCLE  Output 

The  duty  cycle  output  from  DCYCLE  consists  of  printout  only.  Two  formats 
for  pulse  data  are  available  by  option:  a  full  page  printout  (the  same 
as  in  PULSE)  or  a  three  line  printout.  Generally,  most  of  the  pulses 
will  be  printed  with  the  short  format  with  only  one  to  three  full  page 
printouts  per  pulse  sequence.  The  first  line  of  the  short  format  printout  lists 
the:  time  into  the  duty  cycle  of  the  on-signal;  electrical  pulse  on-time; 
off  time  before  and  after  the  pulse;  chamber  wall  temperature  at  start; 
cutoff,  end  of  tail-off  and  at  midpoint  of  on-time.  The  second  line 
lists  the:  total  impulse  factor  to  correct  for  variations  in  energy  losses 
during  a  duty  cycle;  total  fuel  and  oxidizer  flow;  total  impulse;  pressure 
rise  response  times;  pressure  drop  response  time;  and  the  maximum  pressure 
overshoot.  The  third  line  lists  the:  start  and  decay  total  impulse 
bits  interpolated  from  PULSE  tables;  total  impulse  bits  during  pressure 
rise  and  drop  response  periods;  and  mean  pulse  specific  impulse  and  specific 
impulse  efficiency.  Values  are  tabulated  under  column  heading  using 
‘  FORTRAN  coded  names.  For  more  detailed  description  on  these  parameters  ’ 
check  code  descriptions  ir.  Table  4. 

Cumulative  pulse  duty  cycle  performance  data  is  printed  out  on  a  single 
line  at  the  end  of  each  pulse  sequence  and  each  full  page  pulse  printout. 

On  this  line,  cumulative  fuel  and  oxidizer  flows,  cumulative  total  impulse, 
mean  specific  impulse  and  mean  mixture  ratios  of  consumed  propellants 
are  printed. 
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